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TECHNICAL PAPER 


DEVELOPMENT OF IN-SITU STIFFNESS PROPERTIES FOR 
SHUTTLE BOOSTER FILAMENT WOUND CASE 

I. INTRODUCTION 


Recent experiences in formulating and producing appropriate shell stiffness properties for the 
graphite-epoxy wound motor case has necessitated an extension of prediction and verification methods over 
earlier practices. While neither theory nor fiber data was lacking at the inception of the Filament Wound 
Case (FWC) Project, the multitude of compelling structural requirements and the immense geometry of the 
article has rendered some engaging examples for discussion, enticement, and for more comprehensive 
planning on future applications. 

The impetus for substituting filament wound case segments for the proven steel ones was the promise 
of graphite fibers’ high specific strength and stiffness to reduce the vehicle weight such as to increase the 
payload delivered from the Western Test Range by 18 percent. Existing operational constraints were 
imposed and resulting technical challenges identified. The case segmentation was to comply with current 
casting, transportation and erection facilities, but no comparable jointed filament wound pressure vessel 
had been manufactured to provide the necessary confidence to commit the project. Of equal importance 
was the structural requirement that the Shuttle System’s existing load limits were not to be exceeded, but, 
again, no comparable stiffness and dispersion data was available. 

The Solid Rocket Booster Filament Wound Case (SRB FWC) program was preceded by two inde- 
pendent feasibility studies from February through October 1981. Studies concluded that a filament wound 
motor case should be a low technical and schedule risk program. It demonstrated through subscale model 
tests that a segmented FWC was highly feasible and it particularly disclosed the current understanding, 
methodology, and data base for coping with the structural stiffness requirements of the Shuttle systems 
loads. It was immediately perceived that this compliance depended primarily on developing an appropriate 
composite material data base and on developing a brief test program to achieve it. To this end, large bottle 
test specimens were negotiated for developing the A-basis allowables, which included extensional elastic 
constant determination, instead of the standard and more economical 6-in. bottle specimen. Sensitivity 
studies, statistical sampling test, subscale and component verification for math model updating were 
incorporated into the finally negotiated project to minimize major technical problems on the more time 
critical full-scale wound segments. The study could not dispel the suspicion of a scaling effect nor could 
those concerned communicate their concerns - scaling what and why? That was a bridge to come. 

The FWC development has certainly contended with its share of problems, including the in-situ 
membrane response, which almost deluded many participating composite analysts.The extent of this learning 
was demonstrated by the many significant baseline changes of the composite design. What began as a 
material compliance devoid of new concepts and fundamental insights, each new graphite bottle design and 
hydrotest uncovered mistaken computations and stiffness math modeling which persisted during most of 
the development phase. It was not until the third quarter-scale hydro and dynamics tests, when results were 
diligently examined and compared, that the unpredictable responses suggested to the author a new awareness 
of essential phenomena that lead to discovering unexpected parameters which significantly modified the 
data analysis method and modeling techniques. 
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Results presented in this paper are by no means conclusive but are intended to initiate a probe to 
force a solution into this obscure phenomena of filament wound composites which is quite apart from the 
extensive application of graphite cloth. All data and techniques used by the author are presented in the hope 
that others inclined to hack at this elastic code may circumvent any prejudices and perhaps misguided 
intuitions incorporated into preliminary conclusions. Final conclusions will have to wait until the bulk of 
flight-type specimen are manufactured and hydro-proofed. 

At this writing, only one flight-type, full-scale FWC segment has been tested. It represents the largest 
jointed, graphite-epoxy vessel ever built to operate with an internal pressure of 1000 psi. A more striking 
distinction will be achieved when the FWC segments are man-rated. Nevertheless, a number of prescriptions 
for the general development achievements of the elastic compliance are clearly evident. First and foremost 
was the indepth knowledge of composite mechanics, the fabrication experience and mastery of instrumenta- 
tion and meaningful testing demonstrated by the manufacturer through numerous reports, reviews, and 
discussions. Of equal importance is that dynamicists were engaged in the composite description and elastic 
characteristics from the outset of the feasibility study. It follows that as the specific strength of large elastic 
structures increases, the elastic potential energy becomes more critical and the loads and dynamicist role 
more prominent. That role becomes eminent in all phases of design, development and changes when struc- 
tural material is anisotropic. Finally, early Technical Interchange Meetings with associated contractors and 
centers were indispensable for the performance deficiencies revealed that were more suitably assessed and 
resolved by the affected users. 

Special recognition is here due Dr. George McDonough for staffing this effort and for his encourage- 
ment during real time data flow and respective events that led to new basic conclusions that might otherwise 
have been lost to investigators. The author further acknowledges data and reports generated by W. Hill, J. 
York, J. Schutz and C. Kirchner from Hercules, Inc. that were ultimately constructed into Tables 1, 3, 5, 
6, 7, 8, and 9. 


II. FILAMENT WOUND CASE CONFIGURATION 


The solid rocket motor case of the Space Shuttle is segmented to simplify handling, shipping, a-.d 
propellant grain casting. The four filament wound cylindrical segments shown in Figure 1 are field-spliced 
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Figure 1 . Solid rocket motor. 
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by pinning the steel tang and clevis rings as illustrated. Steel-to-steel joints are essentially identical to the 
existing steel motor case and are interchangeable. The composite cylinder ends are joined to the steel rings 
by a similarly pinned tang and clevis arrangement. The length of these segments is defined by field splice 
stations which are common to the steel motor case field joints. The two end domes and external tank 
attach segment are the existing steel configurations and are applicable to either steel or FWC segments. 

< Because the FWC and existing steel motor case must be interchanged on the Solid Rocket Booster and 
must interface with existing Shuttle elements and launch facilities, the end result is that seven steel case 
segments are replaced by four composite segments. 

The composite shell is manufactured from graphite fibers wetted with epoxy and wound on an 
aluminum mandrel in layers. Each layer thickness and fiber direction varies by design to satisfy membrane 
stiffnesses and joint integrity. The membrane construction is dominant and consists of alternating hoop 
and balanced helical laminae. The layup is essentially symmetrical with the mid-thickness which avoids 
coupling between extensional and bending deformation so that the matrix By * 0 forms a more efficient 

structure. The joint ends were thickened with interspersed layers of axially oriented fiber to reduce pin-hole 
bearing stresses and locai deformation. Each case segment is then cured at 390°F through a specified 
schedule. The inner surface is ultimately insulated with a rubber coating and cured at 290°F before casting 
the propellant grain. 


III. COMPOSITE STIFFNESS REQUIREMENTS 


Each FWC segment represents an orthotropic cylindrical shell that was optimized for structural 
strength and stiffness with weight. The shell is essentially subjected to in-plane loading (Fig. 2) so that the 

classical laminate plate theory 111 is a convenient 
basis for defining the symmetrical (or specially 




for extensional and 


Figure 2. Membrane element. 
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for bending deformations. Consequently, all stiffness constraints and requirements imposed on the FWC 
performance were manifested by the Ay and Dy matrices which follow. 
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a) Axial Growth Rate is a critical stiffness constraint during thrust buildup because it excites struc- 
tural components attached forward of the solid rocket booster, the external tank lox dome in particular. 
The thrust buildup rate was fixed so that the axial excursion between the forward and aft skirts was limited 
to 0.6 in. when subjected to combustion pressure of 1004 psig and 3.3SM lb thrust. When the elongation 
resulting from joint tolerances and deformation was deducted, the elongation budgeted to the membrane was 
expressed by 


8 = Le x 


L 

7 





d 


where 


N x 



T 

2»R 


and Nfl - PR 


(3) 


Since the case envelope and loads were constrained by existing Shuttle interfaces, the only opened parame- 
ters were the composite thickness and engineering elastic constants. When these were translated into matrix 
components of equation (1) by the relationships 


1 _ a 22 1 A 1 1 
Ext a 11 a 22" a ?2 E0t A 11 A 22 -A 12 


a 12 _ a 12 
E 6 t E x l A ! l a 22 ~ A^ 2 Xd A 22 


the elongation requirement was appropriately rewritten as 


9. = L e x = L 


N x a 22 ~ N 0 a 12 
a H a 2 2~ a i2 


(4) 


b) Bending Axial stiffness was also specified to avoid the FWC fundamental frequencies from 
coupling with the vehicle control system frequencies of approximately 0.2 Hz. This stiffness was defined 
as the product of laminate variables 

E x t > 10 7 lbs/in. 

or (5) 

t A ll A 22 ' a 12 1 / a 22^ 1q7 • 
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c) Radial Growth of the forward segment was restricted to not exceed the propellant grain tan- 
gential strain capability and was specified by 


R 


AR=Rt0 = - 


N0 

_E 0 



or 


( 6 ) 


N 0 A 1 1 ~ N x A 12 ~ 
.An A22 - Aj 2 . 


<0.58 in. 


d) Structural Integrity of the membrane also implied a stiffness constraint. The tangential unit 
load was much larger than the axial unit load, equations (3), so th»t the hoop fiber ultimate strain was 
related to the hoop fibei ultimate stress, 


Of — Ef €f 
'ult. *L ult. 


or 


(7) 


a f, 


ult. 



N 0 A 1 1 ~ N x a 12~ 
2 

Aj 1 A22 “ Aj 2 _ 


ULT. 


because the hoop fiber strain is identically equal to the laminate tangential strain. 

e) The only requirement imposed on the bending stiffness matrix was that the Dy components be 

defined within the accuracy essential for modal analysis to support dynamics response studies and for 
buckling analysis to sustain water impact loads. 

Having implicitly related the laminate properties to constraints by the above equations, the design 
challenge was to satisfy them all with minimum structural weight through material selections, winding 
geometries and controlled processing. The intermediate graphite fiber was selected because of its high 
specific strength and modulus, and its available data base along with economic considerations. The only 
structural requirement on the epoxy was its ability to bind the fibers together and effectively distribute 
the multi-axial loads. The primary criteria was therefore to select an epoxy that was compatible with manu- 
facturing processes, handling, and operational environments. Performance wise, the lowest epoxy to fiber 
mix ratio contributed the highest specific strength and stiffness in the fiber direction. 

Satisfying the FWC requirements through winding geometries was not so obvious and demands a 
brief review of all winding parameters available to the laminate designer and will be covered in Section IV. 
Even less obvious were the processing controls and independent parameters affecting the lamina which 
makeup the laminate and was the subject of this investigation. 
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IV. LAMINATE DESIGN PARAMETERS 


Adjacent layers having common winding angles are identified as a lamina having orthotropic proper- 
ties with respect to its principal axes (Fig. 3), and are translatable to another set of axes such as the 
cylindrical axes. Because the wound lamina transverse properties in the 2-direction are equal to those in the 
a 3-direction, it is further classified as transversally 

orthotropic and the stress strain relationship, 


°i ICjjl ej 

is characterized by a total of five elastic constants 
which are related to engineering constants as 
follows: 



Figure 3. Wound lamina coordinates. 


C U “ 


Ej -E 2 v \ 2 


C 22 ~ 


e 2 e 1 


E l-E 2 "12 


C 12~ V 12 C 22 


Eo 


C 44 “ C 23 “ 


2 ( 1 + 1 ^ 23 ) 


( 8 ) 


C 55 _C 66' G 12 


C 23 ~ C 22 - : 


-44 


All other matrix components are either symmetrical or zero. 

The FWC is plied in the hoop lamina with the fiber 1 -direction oriented a t a = 90 deg from the 
cylinder x-axis and the balanced helical laminae are oriented at some specified +a and -u angles. Sum of all 
hoop lamina thickness to total laminate thickness, 


H = 



hoop thickness 


Laminate thickness 


(9) 


is a useful design index if all hoop lamina p r operties are identical, otherwise, a modified ratio is necessary. 

To correlate the hoop and helical lamina properties to a common axis, cylinder x-axis, a fourth order 
transformation is used with equations (8) to give the transformed reduced stiffness: 


C H = c n m 4 + 2(C 12 + 2C 66 )m 2 n 2 + C 22 ° 4 

C 12 = (C 11 +C 22 -4C 66 )m 2 n 2 + Ci 2 (m 4 + n 4 ) (10) 
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c 22 = C! ! n 4 + 2(Ci2 + 2C 66 )m^ n‘ + C 22 m 4 

C 66 = < C 1 1 4 C 22 ~ 2C 12 “ 2C 66 )m2 n2 4 C 66 (m4 + ° 4) 

C 16 = < C 1 1 “ C 12 ' 2C 66 )n ’ 3 n 4 < C 12 ~ C 22 4 2C 66^ m " 3 

^26 = ^ C 11 " C 12 ~ 2C 66 ) rnn34 ^ C 12 ' C 22 42 <: 66 > m3n ( 10 ) 

(Cont.) 

^23 = c 12 n ~ 4C 23 m ‘ - 


C 44 = O 44 m 2 + Cjr n 2 


^55 - ^44 4 C 55 m 2 


where m = cosot and n = sina. 

Integrating these hoop and helical plied laminae properties over the laminate thickness in accordance 
with Figure 4, 


N 

Ajj = L 

k=l 


N 

D rt 

k=l 


results into the laminate stiffness components where 
Aji and A 22 are in the cylinder x and 0 axes to 

conform with constraint equations (4). (5), ( 6 ), 
and (7). 

Figure 4. N-plied laminate about mid-plane. Note that the laminate stiffnesses of 

equation ( 11 ) incorporate only three independent 
variables, a, H and t If the designer increases the thickness to satisfy the constraint equations of Section III, 
he will violate the minimum weight requirement. If he increases the wrap angle or the hoop ratio, he satisfies 
the axial growth and bending stiffness but encroaches on the grain and case integrity (Fig. 5). An optimiza- 
tion was indicated during the feasibility study, but until the stiffness property dispersions were developed, 
margins could not be realistically assigned to these constraints. 



( c ij) k U k - z (K-l)l 


(ID 


3 3 
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Figure 5. Stiffness sensitivities. 

A minor infringement on optimum design was that all segments were specified to have a common 
helical angle and hoop ratio so that a single composite system with "A-basis" properties and dispersion need 
be developed. The variation of combustion pressure and external loads acting along the total length of the 
FWC were satisfied by varying segment laminate thickness. The resulting weight increase was not significant 
compared to the benefits of cost and schedules. 


V. LAMINA PRELIMINARY MODEL 


Given that »he fiber, resin and mix have been ideally selected, how might they be modeled to repro- 
duce an experimentally derived lamina elastic response was the next consideration. The Halpin-Tsai [21 
model is extensively used and was favored for this study. The Halpin-Tsai equations not only contain 
parameters that may be adjusted for manufacturing variations, but they reduce the more rigorous model to 
an approximate form consisting of 


M f (l + £ V f ) + M m £ (1 -V f ) 

M ” (1 - Vf) Mf/M m + £ + Vf 

to bound the lower limit properties and the rule of mixtures 


M = M f V f + M m (l -V f ) 


(13) 


to bound the upper limit properties when the constant {-*■»». M symbolizes the elastic property to be 
determined and the known input fiber and matrix (resin) properties are respectively subscripted. Because 
the Halpin-Tsai constant, £, was experimentally derived for each new lamina system, it also compensated 
for ineptly defined constituent properties. This feature was exploited in defining the HBRF-55A resin 
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input properties. Ove^ a dozen neat resin samples were cured, post-cured, and uniaxially tested. The tensile 
modulus spread was between 0.36 and 0.60 Msi and a mean of E m = 0.49 Msi was estimated for the model. 

Other neat resin properties were also estimated because the real strain state was uncertain and the model 
was conveniently forgiving. 

The AS4W-12K graphite fiber properties measured from each lot are listed in Table 1. Orthotropic 
properties of the fiber were estimated from References 4 through 7 for carbon filament (less than 3000°F 
exposure) and noted in Table 2. 

table i . fiber lot properties AS4W - 1 2K graphite 


Lot 

No. 

Modulus MSI 

Pf , 
Ib/in. 

Area 
per tow 
lO^in- 

■■ 

Secant 

Efll 

33.5 

35 



IS 

33 

34 



607-3N 

35 

36.7 

• 1 


607-3P 

34 

35.8 


7.315 

616-4F 

33 

34 

0.0648 

7.054 

61 7-4F 

33 

34.1 

0.0652 

7.066 

6 1 "MG 

33 

34.7 

0.0649 

7.105 

618-3C 

33 

34.7 

0.0647 

7.397 

6 19-3 A 

34 

35.3 

0.0642 

7.530 

62 1-3 A 

34 

34.9 

0.0654 

7.037 

629-3C 

34 

34.5 

0.0652 

7.101 

630-3A 

34 

35.9 

0.0653 

7.049 

632-4B 

33 

35.5 

0.0649 

7.435 

634-4F 

34 

35 

0.0644 

7.547 

634-4G 

34 

35.3 

0.0644 

7.564 

640-4B 

32.6 

34 

0.0646 

7.385 


TABLE 2. ASSUMED CARBON FILAMENT PROPERTIES 


Symbol 

Orthotropic Properties 

Nominal 

E fT 

Transverse Young Modulus, Msi 

2.2 

Gf LT 

Transverse Shear Modulus, Msi 

3.3 

^L 

Longitudinal Poission's Ratio 

0.29 


















To model the series connected constituents, equation (12). a series of simulated lamina tubes were 
tested to determine the Halpin-Tsai constant. Cylinders 4-in. in diameter were wound in hoop only and then 
cured and post cured. Forty-eight specimens were cut to 3-in. lengths and each was fust tested in torsion 
well helow the ultimate and then stressed in tension until failure. Results are summarized in Table 3. 


TABLE 3. LAMINA TUBE TEST RESULTS 13] 


— 
! Symbol 

■ 

Lamina Properties 

— 

Nominal 

1 

l 

Percent C.V. j 

- 

G 12 

Inplane Torsion Modulus. Msi 

0.64 

10 

E 2 

Transverse Modulus. Msi 

0.99 

9 ! 

C2 ULT 

Ultimate Transverse Strain, percent 

0.293 

I 

v t 

Fiber Volume, percent 

51.2 

- 

v d 

Void Volume, percent 

5.9 

— 

1 1 





The resin modulii were degraded for voids. 


E m = K m^- V d> 

(14) 

Gjn ~ E m ’ 2 1 1 ~ v ni ) 

and substituted into the rearranged equation (12) 


Ml(l -V f )M f /M m + V f l -M f 

t= M m- M + v f< M f-Mm> " 5) 

Transverse properties (Table 2) and 4-in. tube ’amina test results (Table 3) were used with equation (15) to 
obtain the Halpin-Tsai constants for the epoxy dominated lamina modulii: 


Young Transverse, = 1,7 
In-Plane Shear. $0 = 2.7 


(16) 


The lamina modulus in the fiber direction and the major Poisson's ratio are parallel connected 
mixtures and comply with the rule of mixtures, equation (13). Correction for void content is unnecessary 
because it is offset by the multi-strain state of the resin. The lamina preliminary stiffness characteristics were 
accordingly modeled by 
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E r E f L V f +E n,"-Vf) 


2 = 0.29 Vj- + 0.35 (1 -V f ) 

E fT <I + l 7 V f ) + E^ 1.7 (1-V,) 
fc2= <>-VflE fT E^tl.7*V, 

Gf LT 0 * 27V f l + G m - 7 < 1 - v f» 
C|: * (l- v f) c r LT . ; C m »--7-‘ V f 

Selected design conditions of 


( 17 ) 


Vf = 0.55 , V d = 0.05 - E^ 33.4 Msi 


were incorporated to predict the FWC lamina elastic properties (Table 4). 

TABLE 4. LAMINA PRELIMINARY PROPERTIES 


Symbol 

Lamina Design Properties 

Nominal 

E l 

Longitudinal Modulus. Msi 

18.6 

E 2 

Transverse Modulus. Msi 

1.0 

G j 

Shear Modulus, Msi 

0.7 

v 12 

Poisson's Ratio 

0.3 


Before winding, the mandrel is covered with a 0.02-in. thick graphite cloth to provide a smooth 
interface and serve as an inner surface protective structure during the insulation coating process. A similar 
surface is provided on quarter scale articles by a 0.02-in. thick glass cloth to scale down its contribution to 
the total stiffness. In each case, the cloth is assumed to act as a lamina 0.01-in. thick ai the woof direction 
and 0.01 -in. lamina in the warp direction. Respective laminae properties are listed in Table 5. 
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TABLE 5. CLOTH LAMINA PROPERTIES 


Symbol 

Elastic Properties 

Glass 

Graphite 

E 1 

Longitudinal Modulus. Msi 

6.6 

17.4 

i E-? 

Transverse Modulus. Msi 

2.0 

1.6 

G | 

Shear Modulus. Msi 

0.76 

071 

Vp 

: 

Poisson’s Ratio 

0.28 

■ — J 

0.23 


VI. LAMINATE PERFORMANCE TESTS 


As is frequently the case, developing experimental evidence to verify design predictions is the most 
technically exacting and programmatically demanding phase of many high performance projects. Testing 
tasks are further compounded when a statistical basis is imposed on subscale specimen. 

In reviewing constraint equations (5) through (7). it was noted that all but one could be directly 
verified by measuring the axial and hoop strains during hydrotest. It was then believed that since the most 
common procedure for establishing the A-basis strength of a laminate system was to hydroburst numerous 
6-in. bottles, appropriate instrumentation might ride on larger bottles to back out. statistically [8) the most 
significant elastic constants using acceptance tested fiber of 1 0 or more lots. 

The specimen selected was a 20-in. diameter bottle based on expediency and cost. The bottle was 
designed to be wound on a sand mandrel covered with 0.06 in. cured silicone rubber. End closures were 
wrapped with the baselined 29-deg helical angle fibers, though these end domes were intended for a 25-deg 
wrap angle. The winding sequence was helical-hoop-helical-hoop-helical for a total thickness of 0.181 in. 
and stress ratio of helical to hoop of 0.29. Layer thicknesses and properties varied according to the fiber 
lot used. Bottles were cured and post-cured in a 6 to I nitrogen-air environment. The bottles were instru- 
mented with 3 linear variable differential transformers (LVDT) to measure axial displacement and three 
linear potentiometers (LP) to measure hoop growth. Strain gauges were also applied along the principal axes 
for comparison. 

Manufacturing complication with the compromised mandrel emerged at once. Though the helical 
angle was only 4 deg in excess of the available 25-deg mandrel, fibers began to siip at the dome ends until 
the winding operator learned to manually resist the slipping wetted tow. This operator influence was not 
expected to restore the programmed uniform winding tension on the fiber. The second objection of the 
sand mandrel was that the expansion of the full scale aluminum mandrel during cure was not being simu- 
lated. Both objections were believed to cause slack fibers resulting in low fiber strength and large dispersions. 
Nevertheless, this series designated as ABA specimen, provided an excellent opportunity to assess instru- 
mentation influences and laminate response verification limitations. 

On instrumentation, LVDT proved to be far superior for measuring axial strains consistently and 
accurately than strain gauges. However, measuring LVDT’s was very sensitive to membrane discontinuity 
stress waves which were noted to propagate over 30 thicknesses from the dome-cylinder intersection. The 
LP’s consistently measure hoop strains about 4 percent low. 
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Strength and strain dispersions were exorbitantly large though the last 34 bottles, identified as ABA2 
series, were an improvement over the first batch because the operator influence improved. Properties and 
hydrotest data for these ABA2 specimen are documented in Table 6. Later they were abandoned as “A- 
basis” data, but are included here for scaling trends. 

In the meantime, the quarter-scale model case segments were being developed as deliverable articles 
to verify the 1 WC dynamics response for on-pad. liftoff, nominal flight, and burnout conditions. The seg- 
mented case was 350 in. long with an inside diameter of 36.25 in. measured at the membrane. Metal and 
composite joints simulated those of the full-scale article and the helical wrap angle was that of the pre- 
liminary baseline. ±29 deg. 

A number of these segments were scheduled to be hydroproofed to verify the membrane response 
as well as joint integrity. This test series was identified as DV-36 specimen. To insure that the composite 
experienced the same thermal exposure as the full scale, a post-cure was imposed on the segment after it 
was removed from the mandrel to simulate liner vulcanizing environment. It was here discovered that this 
post-cure caused intolerable radial dimensional changes attributed to stress relief. All subsequent segments 
were specified to be stress relieved before machining composite ends for joint adapter rings to avoid 
locked-in discontinuity stresses during liner cure. 

Based on the first three hydroburst specimens, it was noted that the average fiber strength increased 
10 percent over the ABA2 specimen and was concluded that the DV-36 series was more representative of the 
manufacturing and cure processes because of the improved mandrel and the increased lamina layers. The 
“A-basis” allowable test emphasis was switched from the 34 specimen ABA2 series to the 20 specimen 
DV-36 series. Manufacturing and hydrotest data for 18 specimens are listed in Table 7. 

The quarter-scale dynamics test articles were manufactured and hydroproofed before delivery and 
results are listed in Table 8. Though these eight segments belong to the DV-36 family of specimens, they 
were characterized separately in this study for comparison. 

The laminate baseline was updated at the Critical Design Review. Quarter-scale segments of the 
deliverable structural test article (STA) were manufactured with the updated laminate system. A variety of 
structural tests were performed on this quarter-scale series designated as STA Jr. and results from hydrotests 
are noted in Table 9. 

Only three ful-scale test articles have been manufactured and hydrotested. These articles represent 
the aft segment of the FWC and only the DA003 is of the flight-type composite system. Data reported in 
Table 10 was considered preliminary at the time of this investigation. 

In all the above tables, the fiber volume fraction was calculated from 


(1-V d )(l-W m ) 
f '(l-W ni ) + W m p f /p m 


( 18 ) 


where the void fraction was assumed. V d = 0.05. and the resin ratio by weight was controlled at W m = 
0.327. 
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TABLE 6. ABA2 PROPERTIES AND HYDROTEST DATA 






Thicknesses, inch 

— 

Percent Strains 
Measured 

Sp. 

No. 

Fiber 
Lot No. 

v f 

Percent 


Helical 

Layer 

Hoop 

Layer 


Load. N x 
lb/in. 

V d 

Percent 

Total 

e * 

e 0 * 

1 

616-4F 

55.8 

5.6 

0.426 

0.0245 

0.1768 

6279.5 

0.1411 

0.9838 


1 

54.4 

6.5 

0.0446 

0.0257 

0.1852 

6276.9 

0.1241 

0.8469 

3 


54.0 

6.1 

0.0446 

0.0257 

0.1852 

6276.9 

0.1291 

0.8635 

11 

617-4G 

53.1 

5.1 

0.0429 

0.0247 

0.1781 

6279.1 

0.1369 

0.9595 

12 

6174G 

56.2 

4.6 

0.0429 

0.0247 

0.1781 

6279.1 

0.1357 

0.9882 

13 

617-4G 

53.8 

5.1 

0.0429 

0.0247 

0.1781 

6279.1 

0.1340 

0.9772 

14 

62 1-3 A 

52.8 

5.7 

0.0426 

0.0245 

0.1768 

6279.5 

0.1250 

0.9540 

15 

62 1-3 A 

56.3 

5.6 

0.0426 

0.0245 

0.1 768 

6279.5 

0.1475 

0.9750 

16 

62 1-3 A 

51.9 

6.2 

0.0426 

0.0245 

0.1768 

6279.5 

0.1210 

0.9374 

17 

6174F 

52.5 

4.2 

0.0428 

0.0246 

0.1776 

6279.2 

0.1364 

1.0286 

18 

6174F 

53.7 

5.1 

0.0428 

0.0246 

0.1776 

6279.2 

0.1176 

0.9419 

19 

607-3P 

53.5 


0.0440 

0.0253 

0.1826 

6277.7 

0.1335 

0.9143 

20 

618-3C 

53.4 


0.0446 

0.0257 

0.1852 

6276.9 

0.1319 

0.8767 

21 

618-3C 

56.0 


0.0446 

0.0257 

0.1852 

6276.9 

0.1182 

0.9209 


607-3N 

53.9 


0.0455 

0.0262 

0.1889 

6275.7 

0.J322 

0.7947 

23 

607-3N 

52.0 


0.0455 

0.0262 

0.1889 

6275.7 

0.1192 

0.7847 

24 

6174F 

54.0 

4J 

0.0428 

0.0246 

0.1776 

6279.2 

0.1368 

0.8667 

25 

607-3P 

53.3 

5.7 

0.0440 

0.0253 

0.1826 

6277.7 

0.1268 

0.8893 

26 

607-3P 

52.2 

9.5 

0.0440 

0.0253 

0.1826 

6277.7 

0.1268 

0.8775 

-»? 

618-3C 

52.8 

7.8 

0.0446 

0.0257 

0.1852 

6276.9 

0.1304 

0.8480 

28 

618-3C 

53.6 

6.1 

0.0446 

0.0257 

0.1852 

6276.9 

0.1293 

0.8500 

29 

607-3N 

54.1 

5.8 

0.0455 

0.0262 

0.1889 

6275.7 

0.1226 

0.8499 


607-3N 

53.2 

6.9 

0.0455 

0.0262 

0.1889 

6275.7 

0.1135 

0.85S5 

31 

629-3C 

56.4 

5.9 

0.0430 

0.0247 

0.1784 

6279.0 

0.1158 

0.8400 

33 

629-3C 

57.4 

4.2 

0.0430 

0.0247 

0.1784 

6279.0 

0.1364 

0.8864 

34 

629-3C 

59.2 

4.9 

0.0430 

0.0247 

0.1784 

6279.0 

0.1886 

0.8201 

Mean 

54.2 

5.98 

0.04374 

0.02517 

0.1815 

6277.3 

0.1318 

0.8973 

Percent C.V. 

3.2 

19.6 

2.4 

2.4 

2.4 

0.044 

10.7 

6.3 


3 Layers at +29 deg balanced 

3 Layers at - 29 deg Helical 

2 Layers at 89.4 deg. Hoop 
p m = 0.0437 lb/in 3 

*C o r :ected to strain gauge 


Hydro press. P = 1 267 psig 
H = 0.2773 avg. 

E f =33.6 MSI, C.V. = 2.1% 
l.D. = 20 in. 
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TABLE 7. DV-36 “A-BASIS” PROPERTIES AND HYDROTEST DATA 


SP 

No. 

Fiber 
Lot No. 

V f 

Percent 

Specimen Thickness, in. 

Load N x 
lb/in. 

Percent Strains 
Measured 



Total* 

Calc. 

Total 

Meas. 

Per L 
Helical 

ayer 

Hoop 

e x 

e d 

1 

— 

607-31 

54.9 

' 

0.0491 

0.0284 

0.2816 

0.288 

9980.8 

0.9577 

0.2001 


607-31 

54.9 

0.0491 

0.0284 

0.2816 

0.288 

9980.8 

0 8756 

0.2106 

4 

6 19-3 A 

55.8 

0.0505 

0.0293 

0.2899 

0.291 

9978.5 

0.9244 

0.1868 

5 

617-4G 

55.5 

0.0477 

0.0277 

0.2739 

0.282 

9982.9 

1.0169 

0.1916 

6 

630-3.4 

55.4 

0.0472 

0.0273 

0.2707 

0.273 

9983.8 

1.0033 

0.2033 

7 

629-3C 

55.4 

0.0477 

0.0276 

0.2736 

0.278 

9983.0 

0.9777 

0.1894 

8 

6174G 

55.5 

0.0477 

0.0277 

0.2739 

0.2’4 

9982.9 

1.0166 

0.1958 

9 

629-3C 

55.1 

0.0476 

0.0276 

0.2732 

0.284 

9983.1 

0.9776 

0.2020 

11 

6 19-3 A 

55.4 

0.0505 

0.0293 

0.2899 


9978.5 

0.9094 

0.1861 

12 

632-4B 

55.2 

0.0499 

0.0289 

0.2803 

0.269 

9979.5 

0.9472 

0.1869 

13 

6344F 

55.3 

0.0506 

0.0293 

0.2903 


9978.4 

0.9431 

0.1874 

14 

6344F 

55.3 

0.0506 

00293 

0.2903 


9978.4 

1.0160 

0.1957 

15 

629-3C 

55.1 

9.0477 

0.0276 

0.2736 


9983.0 

0.9902 

0.1882 

16 

6344G 

55.4 

0.0507 

0.0294 

0.2910 


9978.2 

0.9601 

0.1979 

17 

6344G 

55.4 

0.0507 

0.0294 

0.2910 


9978.2 

0.9603 

0.1941 

18 

630-3A 

55.0 

0.0472 

0.0273 

0.2707 


9983.8 

0.9689 

0.1898 

19 

630-3A 

55.0 

0.0472 

0.0273 

0.2707 


9983.8 

0.9928 

0.2023 

20 

629-3C 

55.1 

0.0476 

0.02 7o 

0.2732 


9983.1 

0.9725 

0.1984 

Mean 

55.3 

0.04885 

002830 

0.2802 

0.281 

9981.2 

0.9672 

0.1948 

Percent C.V. 

0.42 

2.90 

2.94 



2.91 

2.6 

0.02 

3.8 

3.5 


•Thickness excludes .02 in. glass cloth 


p m - 0.0437 for 1 and 2 
P m = 0.0444 for 4 through 8 
P m = 0.0437 for 9 through 20 
3 layers at +29 deg ) balanced 
3 layers at - 29 deg ^ helical 
2 layers at 90 deg hoop. 

1 layer glass cloth 0.02 in. thick 
Mean E^ = 33.78, 1.12 percent C.V. 

Hydro press, p = 1110 psig 


I D. = 36.25 in. 

Calculated effective thickness: 
t e = t + tg Eg/Ef = t + 0.0037 in 
H = 0.3029 
Void, Vj = 5 percent 
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TABLE 8. Q. S. DYNAMICS ARTICLES PROPER'i 'ES 


■ 

Fiber 

Lot Number 

V f 

(%) 

■ 

Measured 

Thickness 

(in.) 

Percent Strains 
Measured 

e x 



e 6 

1-F 

1 

61.9 

0 

■ . 

mm 

0.8300 

1-FC 

1 1 

61.3 

1.0 

1 

W&SiM 

0.8280 

1-AC 

607-31 

60.5 

0.5 


0.1877 

0.8900 

1-A 

607-31 

56.6 

1.6 


0.1927 

0.8200 


607-3N 

58.0 

1.7 

0.306 

0.1781 

0.8500 


617-4F 

55.8 

1.6 

0.285 


0.9440 


7-4F 

62.5 

0.5 

0.279 


0.9430 


’-4F 

60.2 

0.8 

0.280 

0.1458 

1.0300 

Mean 

56.9 

0.96 

0.294 

0.1846 

0.8919 

Percent C.V. 

4.1 

61.0 

3.8 

9.8 

7.8 



Helical 

Hoop 

Graphite thickness 0.2816 in. 

Angle, deg 

±29 

90 

Glass cloth 

0.02 in. 

No. Layers 

8 

3 

Total thickness 

0.3016 in. 

Layer Thickness 

0.02455 

0.0284 

I.D. = 36.25 in. 



Unit load, N = 8991 lb/in. H = <*.3025 

Hydrotest pressure, p = 1000 psig 

TABLE 9. STA JR. PROPERTIES AND HYDROTEST DATA 


SP No. 
DV-36 

Measured 

Thickness 

Hydro 

Pressure 

psig 

Load 

N x 

lb/in. 

Percent Strains 
Measured 

e x 

e 6 

-25 


1102 

9905.5 

0.0650 

0.9554 

-26 

0.335 

1108 

9959.5 

0.0400 

0.9560 

-27 


1108 

9959.5 

0.0424 

0.9904 

l 

00 


1114 

10013.4 

0.0441 

1.0086 

-31 


1114 

10013.4 

0.0442 

0.9981 

-33 


1110 

9977.4 

0.05612 

0.9443 

-35 


1114 

10013.4 

0.04322 

1.0010 

-36 


1118 

10049.3 

0.04838 

0.9753 

-38 


1110 

9977.4 

0.05548 

0.9186 

Mean 

1110.9 

9985.4 

0.04882 

0.97197 

Percent C.V. 

4.0 

4.0 

16.0 

2.9 



Helical 

Hoop 

Glass cloth 0.02 in. thick 

Fiber Lot No. 

640-4 B 

630-3A 

Voids, Vj = 5 Percent 

Angle, deg 

±33.5 

90 

Total thickness, t = 0.314 in. 

No. Layers 

12 

5 

H = 0.245 

Layer thickness 
Fiber Volume, Vf 

0.0185 

55.0% 

0.0144 

55.3% 

I.D. = 36.25 in. 
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TABLE 10. FULL-SCALE SEGMENT PROPERTIES AND TEST DATA 




Thickness* 

Load 

Percent Strains 


v f 

Percent 

inc 

les 

N x 

Measured 

Sp. Number 

Calculated 

Measured 

lb/in. 

e x 

e 0 

DA001 

55.2 

1.181 


32453 

0.05144 

0.9670 

DA002 

55.1 

1.143 


33182 

0.05750 

0.9850 

DA003 

58.0 

1.308 

1.447 

37962 

0.05043 

1.0086 


Winding Lay up. 

DA001 

DA002 

DA003 

Hoop angle, deg 

90 

90 

90 

No. layers 

4 

4 

4 

Thickness 

0.0312 

0.0306 

0.03048 

No. layers 

9 

9 

11 

Thickness 

0.0156 

0.0153 

0.01524 

Helical angle, deg 

±31.85 

±31.85 

±33.55 

No. layers 

28 

28 

6 

Thickness 

0.03195 

0.0313 

0.0156 

No. layers 



14 

Thickness 



0.0312 

No. layers 



6 

Thickness 



0.052 

No. layers 



2 

Thickness 



0.0312 

Hoop ratio, H 

0.2283 

0.2288 

0.234 

Fiber Lot No. 

634-4F 

636-4G 



-4G 

-4L 


Hydro pressure, psig 

901 

921 

1051 


♦Includes 0.02 in. graphite cloth 
I.D. = 145.23 in. 
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The axial unit loading experienced by the surface mounted gauges was approximated by 


n 

* 2R + 1 


(i9) 


and the thickness of the glass cloth used in quarter-scale specimen was reduced by its stiffness contribution 
to provide an effective laminate thickness of 


tM'-MjEj/Ef 


where t' is the laminate thickness minus the glass cloth thickness and subscript g refers to glass cloth proper- 
ties. 

As previously mentioned, hydrotests will verify most of the project compliance except the bending- 
axial stiffness and the diagonal components of Ay and Dy matrices. Equation (5) may be directly verified 

by uniaxial tests of dog-bone specimens from full-scale articles. One set of such tests has been performed on 
DA003 membrane. 


E x • t = 7.85 x 1 0 6 Ibs/in. 


More tests are scheduled on this flight-type membrane which will include uniaxial tests in the hoop 
direction and beam tests on specimen with hoop and axial orientations. Results of these tests will, of course, 
be used to correlate hydrotest analyses. 


VII. LAMINATE VERIFICATION ANALYSIS 


If the laminate performance of a flight-type segment has been verified directly by hydrotest and 
was accepted, why labor over updating the math model? The primary purpose for developing a fine-tuned 
elastic model was to predict the response for “all” FWC operational loading condit’ons such as static firing, 
on-pad cantilevered mode, liftoff and twang transients, ascent flight and reentry through splash-down 
environments. Included among the loading conditions imposed on the FWC is the acceptance hydroproof 
test which is comparable to foregoing development hydrotests. If the math model cannot reasonably predict 
principal strains for a hydroproof test, then there is serious concern that the manned FWC will comply with 
critical flight environments. It would seem that the bonanza of composite development data represented by 
Tables 6 through 10 should provide an adequate basis for formulating required laminate math models. 

The laminate preliminary stiffness components were forthrightly calculated with estimated average 
lamina properties, Table 4 and equation (11). Results for the six different development specimens are noted 
in Table 1 1. These stiffness components cannot be verified directly but may be backed out from a variety 
of tests relating stresses and strains. Instrumenting all scheduled hydroproof articles was not only a prudent 
choice but assured maximum data return. Hydrotests were performed with closed ends so that unit load 
equations were reduced tc 
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N 


x 


PR 

2 


Nfl 

2 


and because both laminate layup and shell geometry are symmetrical, laminate shear is nulled which further 
reduced constitutive relations of equations (1) to two biaxial strain equations, 


N x~ A ll e x + a 12 € 0 


2N x“ a 12 € x + a 22 e 0 


( 20 ) 


involving three elastic properties. Here was the analytical challenge of the hydrotest program, how to 
decrease the infinite possible sets of equation (20) solutions to the most likely one. 

TABLE 11. ESTIMATED LAMINATE PROPERTIES 
(Lamina Modeled from Section V) 


Parameters 

QS1,2 

DV-36 

ABA2 

DA001,2 

STA. Jr 

DA003 

Aj j Mppi 

2.453 

2.439 

1.568 

9.745 

2.347 

10.300 

A 2 2 Mppi 

2.119 

2.111 

1.230 

7.452 

2.087 

8.549 

Aj2 Mppi 

0.672 

0.667 

0.438 

3.253 

0.866 

3.898 

a 66 MPP‘ 

0.787 

0.783 

0.510 

3.716 

0.987 

4.367 

Di j Kpi 

19.97 

19.66 

5.078 

1130.0 

19.69 

1468.0 

D 22 Kpi 

11.75 

11.59 

2.121 

834.0 

14.28 

1253.0 

Dj 2 Kpi 

5.29 

5.21 

1.415 

365.0 

6.96 

530.7 

d 66 K P J 

6.12 

6.02 

1.612 

417.0 

7.90 

603.9 

81 

5.463 

5.419 

7.621 

17.164 

10.783 

19.639 

g 2 Mppi 

1.121 

1.117 

0.640 

3.821 

1.084 

4.372 


To start, equations (20) were rewritten to separate laminate stiffness parameter from response 
parameters, predicted or obtained from hydrotest, 


eg 2 A n -A 12 
e x A 22 “ 2A 12 

2 

N x _ A 11 A 22~ A 12 
60 2Aj i - Aj 2 


( 21 ) 
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The second expression might just as easily be related to the hoop unit load. It turned out that these 
two ratios conveniently and completely characterized the behavior of a pressurized composite membrane 
and are here coined as g n -charaeteristics or constraints. Calculated stiffness properties, Table 1 1 , were used 

to predict the g n -characteristics of equation (21) to compare with those determined from hydrotest data. 

T ... approach was adequate for comparing between predictions with an individual test. To derive these 
characteristics from multiple data of a particular configuration, a curve-fitting method was used with Tables 
6 through 9. 

It was noted that both expressions in equation (21) are of the form 


0 = y i/ x i 


and the problem was to calculate the “best fit” y^ and x^ from n-paired observations Using least square 
method, (3 is chosen to minimize 


I 


i=l 


tyj-0Xj] 2 


Differentiating with respect to 0 and setting to zero yields the normal equation 



or 



and 



( 22 ) 


Results from hydrotest specimen are listed in Table 12. 

Since hydrotest provides the g-characteristics of equations (21). a .er equation is required to 
resolve the three stiffness properties of the test article. That equation wa* fon.,ulated [8] from the principle 
of lease square where the sum of the square of the differences between the observed, Ay, and the estimated, 

Ay, parameters results is the smallest possible. This principle is expressed in functional form by 

f(Ajj) = (A U - Aj t ) 2 + (A 22 - A 22 ) 2 + (A 12 - A 12 ) 2 = min. (23) 

The estimated parameters were derived from preliminary equations (17) whose values are listed in Table 1 1 . 
20 
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TABLE 12. HYDROTEST DERIVED CHARACTERISTICS 
(Tables 6 through 10) 



QS 1,2 

DV-36 

ABA2 

D001,2 

STA. Jr. 

DA003 

No. Bottles 

8 

18 

26 

2 

9 

1 

81 

4.77 

4.96 

6.77 

17.91 

19.35 

20.0 

g 2 . M PPi 

1.00 

1.03 

0.696 

3.3o3 

1.026 

3.7b4 


Using the Lagrange multipliers, X n , the functions' equation is minimized while satisfying the con- 
straints, g n , by the expression 


3Ajj dAy tf (A ij ) 2Z X n 0 ■ (24) 


Rewriting the constraint equations (21), 


01 - 2Aj j -gj A 2 2 + Aj2 (2gj - 1 ) = 0 

2 

0 2 = 2g 2 A u - A]] A 22 -£2 a 12 + a 12 =0 


(25) 


and substituting with equations (23) into equation (24) gives 


-[(Ajj ~ Aj j)“ + (A 22 ~ A 22 )~ + (A] 2 ” Aj 2 )^ + Xj(2Aj j " 8lA 22 + a j 2 (2gj - 1)) 
d 

+ X 2 (2g2 Aj j - Aj j A 22 “ g2 a i 2 + A^ 2 )] = 0 • - ) 


Making equation (26) stationary with respect to the observed variables. Ay, 


- -2(A n -A n ) + 2Xi + 2X 2 g 2 ~X 2 A 22 - 0 

^22 = " 2 ^ 22 ~ A22) * Xl?I ' ^ Al 1 = ° 

= -2(Aj 2 - Aj 2 ) + X 1 (2gj - 1 ) - X 2 g 2 + 2X 2 Aj 2 = 0 
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and solving for \| and X 2 from the first two and substituting into the third gives the desired third equation 
(1 - -gl> Aj ! + [Aj 1 (2gj - 1) - gj g2 - 2A 12 J Aj j + (2gj + 2) A n A J2 


- I2gj Aj j + 4A 22 + 2gj g 2 l Aj 2 + (4 +■ gj) Aj 2 A 22 + [A 22 (1 - 2gj) 


(27) 


-4gj 83 ”®1 A 12^ A 22 + (-Si ~ 1) A 2 2 — gj ®2 l A l 1 + 4A 22 + “ A 12^ 


The Newton iterative method (9] was used to solve nonlinear algebraic equations simultaneously by 
expanding in Taylor series about a first approximation and keeping only the first, such as, 


A .. = A 0) +6 U) A ? = a ? (l) + V -5*. l > 
a ij a ij °u » A ij a ij - a y °ij 


etc. 


Substituting into equations (25) and (27) and solving for 6:: and then adding these values to the first 

J (i) 

approximation to yield a second approximation, the procedure is iterated until 5^ converges. The three 
simultaneous equations written in final form are 


25 j j - gj S 22 + (2gj - 1) 6 J2 — 2a j j + gj a 22 + (1 - 2g j ) aj 2 


(2d) 


(2gi -a-> 2 ) 5j j — a j j 5 22 + (2aj 2 -gT)6j 2 = ajj Aii- 2g-> ajj +g-»aj-»-a| , 2 


(29) 


l(2g, - 1 ) (A, , - 2a, , ) - g, g 2 -2A 12 + 2a 12 (g, + 1)] 5, j + ((2g, - 1) (2a 22 - A 22 ) 


-4gj g 2 -gj Aj 2 + aj 2 (gj +4)) 5 22 + t— a i 1 (gj + 1) - 2gj Ajj -4A 22 ~2gj g 2 
+ a 22 (4 + gj)l 5 j 2 = ajj (2g- l)-ajj (Ajj (2g- l)-gj g 2 -2Aj 2 j 


~ a ll a 22 (2g l + 2) + aj 2 (2gj Ajj + 4A 22 + 2gj g 2 )-a, 2 a 22 (4 + gj) 


-a 22 IA 22 (1 — 2g j ) — 4g | g 2 -gj A, 2 J -a 22 (2gj - 1 ) - gj g 2 Ajj 


4 *1 g 2 A 22 ” 2g I g 2 A 12 


( 30 ) 


Equations (28) through (30) are of the form 
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P| **11 + °-l S'*' - * + r l fip-S) _ 

p-» 6j 1 + q-« 5 -»t + r-iip = s-i 

P3 5 1 1 + 93 5 22 + r 3 6 l2 = s 3 

where 6 jj was solved by Cramer’s rule through the following expressions 
Uj = Pj (Q2 r 3 ~ r 2 93 ) ” P2 ^t r 3" r l < 13> + P3( < ll r 2 -r l 92* 

U 2 = sj (q2 r 3 ~ r 2<l3>-S2 <9j r 3~ r l 93> + s 3 (9i r 2 -r l 92> 

(31) 

U 3 = P 1 (s 2 r 3 ~ r 2 s 3^ " P2 ^ s 2 r 3 " r l s 3* + P 3 < S 1 r 2 " r l s 2* 


U 4 = Pj (qT S 3 - St q 3 > - Pt (42 r 3 _ S 1 <* 3 ) + P 3 (9l s 2 ~ S 1 92* 


5j I - Ut/Uj , 5n-U^/U| , - U^/Uj 


(32) 


In using the above outlined approach, the “most probable” observe*.' properties for each set of 
hydrotest articles may be determined if the “estimated” and “observed" are >ust a few percentages off, 
otherwise, th _• estimated modeling should be open to doubt. Table 13 lists the calculated “observed” 
properties for all six sets of hydrotests and Tab;e 14 provides the comparison in percentage between the 
estimated and observed, which at a glance was bewildering. Here were over a hundred articles built and tested 
to the same specifications, but the prediction disparity ranged 40 percent under to 30 percent over. Further- 
more. no pattern was discemable as to article size or response dependency. Since the laminate is defined 
solely by the layup geometries of laminae layers, differences between observed and estimated must be com- 
pensated by adjusting the lamina elastic models. That was like dumping another eight unknown variables 
with only two characteristic equations to resolve them: the observed laminate properties of Table 13 were 
suspect by the resulting gross disparity: therefore were not counted. 

TABLE 13. OBSERVED LAMINATE PROPERTIES 
[Equations (25) and (27)1 


Laminate 

Stiffness 

QS 1,2 

DV-36 

ABA2 

D001.2 

STS. Jr. 

DA003 

An Mppi 

2.425 

2.420 

1.570 

9.723 

2.349 

10.273 

Att Mppi 

1.897 

1.951 

1.323 

6.568 

2.005 

7.365 

Aj 2 Mppi 

0.492 

0.542 

0.464 

2.820 

0.904 

3.250 
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TABLE 14. DEPARTURE OF ESTIMATED FROM OBSERVED PROPERTIES 

(Percent) 



QS-1,2 

DV-36 

ABA2 

D001.2 

STA. Jr. 

DA003 

A ll 

1.1 

0.7 

0 

0.2 

- 0.1 

0.3 

A->-> 

11.7 

8.2 

- 7.0 

13.5 

4.1 

16.1 

Ap 

36.6 

23.2 

- 5.6 

IS , 

- 4.2 

1S.4 

«1 

14.4 

9.5 

12.5 

-4.1 

-44.1 

-1.8 

§2 

12.1 

8.4 

- 7.5 

13.6 

5.6 

162 


Surely not all lamina properties were causing gj and g-« characteristic disparity and with some well- 

grounded judgement, unknown variables may be reduced to a manageable few. To gain this insight, a 
sensitivity analysis of the g-characteristics with respect to lamina elastic properties was developed. Appendix 
A. 


The sensitivities of the six sets of specimen were noted to vary slightly with geometry but the 
following generalizations are applicable to all. 

The g-» characteristic is the most difficult to adjust. It is most sensitive to the hoop lamina longi- 
tudinal modulus and moderately sensitive (an order of magnitude) to the lamina transverse and inplane 
shear modulii. It is least sensitive to the helical lamina longitudinal moduiur and hoop lamina transverse 
modulus. 


The g| characteristic is most sensitive to longitudinal modulii of both helical and hoop lamina. It 

is one order of magnitude less sensitive to the helical transverse modulus and least sensitive to the hoop 
lamina transverse modulus. 


by 


It appeared that g-» must be adjusted first by modifying the hoop lamina longitudinal modulus (Ej) 


AS g-. 

18.6 I = — 

100 x l8.6(3g 2 /dE,) 


(33) 


because it is less sensitive to adjustments of g| which follows. This was performed for each set of specimen. 
What resulted was that the g 2 adjustment intensely aggravated the disparity of gj characteristic. The new gj 
was adjusted similarly to equation (33). Recalling that g] is most sensitive to the longitudinal modulus of 
the helical lamina, it was modified as noted in Table IS. 


In spite of the directness of this approach, a number of serious difficulties turned up that were 
contrary to experience. While some modification of the in-situ lamina modulii was expected, a 50 percent 
reduction of the fiber dominated modulus is not conceivable, especially in the hoop lamina - the mechanism 
is just not there. That one alone was sufficient cause to dismiss this mode of modification and other 
difficulties became academic. 
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TABLE 1 5. LAMINA TRIAL COMPLIANCE WITH TEST g, AND g 2 


Modified Parameters 

K3B 

DV 

ABA2 

DAI. 2 

STA. Jr. 

DA003 

g-> compliance 







Change Ej, msi 

16 

17 

20 

15 

17 

15 

Result gj percent 

60 

53 

9 

765 


-374 

gj compliance 







Change Ej, percent 

16 

15 

*> 

61 

I- 9 

- 15 

Ej, msi 

15.5 

15.8 

18.2 

7 26 

wmsm L.' 

21.5 


No doubt the lamina properties must be modified for in-situ conditions, but results must pass the 
test of experience however limited at this time. 


VIII. IN-SITU LAMINA MODELS 


In the manufacturing and curing processes, laminate layers experience a wide gradation of polymer 
changes and imposed contact strains (winding tension, cure and thermal shrinkage, etc.). Intuitively, one 
might expect the resin connected properties to be especially affected with hoop layers unlike helical ones 
and thick layers unlike thin ones. Alerted by this possibility and guided by results of Section VII and insights 
from the sensitivity analysis, the approach for modifying the lamina model to satisfy test data was obvious, 
first reduce the resin connected modulii as related to manufacturing phenomena and then make up tire 
remaining disparity with fiber connected modulii adjustments. 

The first step was crucial to the success of the model and it triggered many collected thoughts and 
observations of tow winding and resin distribution. In the beginning, the prepared resin is checked for con- 
sistency but its tackiness changes with time. The stiffer resin may affect the fiber coating and change the 
resin-to-fiber ratio and it may also resist compaction and coalescences of adjacent tows. Consequences would 
be to decrease the lamina modulii and increase the article thickness over that predicted. This occurrence, 
though ignored, is well known to the industry. But in this investigation it provided a baas for scaling - 
TIME scaling of the winding process related to resin content and bonding. 

Time parameter rightfully belongs in the polymer chemistry domain, nevertheless, it may be roughly 
related to some available geometric properties of the different set of articles such as 


E 2 ,G 12 - f(N,t.S) 


(34) 


where N, t are number of layers and thickness of a laminate and S refers to the scale size of the article. 
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Another phenomena affecting the resin con- 
nected moduiii, not time dependent, was realized 
here. Winding tension on the tow (Fig. 6) is fairly 
constant for all size articles, wrap angles, and layers. 
Yet the surface geometry can be shown to vary 
which has a distinct effect on the tow spread. Assum- 
ing the general case of helical winding, the mandrel 
profile near the tow tangent point may be approxi- 
mated as an ellipse with R and *‘b” minor and major 
semi-axes respectively. The equilibrium equation 
normal to the tow is 

d0 dd 

N r r d0 - (F + dF) sin — — F sin — = 0 


or 

N' r d0 - Td0 = 0 

and N r = F/r (35) 


FREE 



Figure 6. Winding tension effect on tow. 


The physical significance of this expression is that for constant tow tension, the normal unit load 
acting on the tow is inversely proportional to the radius of curvative “r.” Since the ability of the tow to 
spread (flatten) and touch the adjacent tow* is related to the normal load, then the extent of coalescences 
and voids (gaps) between them is dependent on the radius of curvature at the tow lead-in tangent point 
(Fig. 6). 

The elliptical profile may be better expressed in winding terms such that the major semi-axis is 
b = R/tan a where a = helical wind angle. For a hoop wound lamina, b = R so that tan a = 1 and the radius 
of curvative is identically the niadrel radius R and the normal unit load is inversely proportional to the 
mandrel radius by equation (35 ). For the helical wound lamina, the radius of curvature is 


r = (x - tan 4 a + y-J 


3/2 


R 


1 “» 

** tan* a 


for a < 45 deg. At the point of interest, x = 0 and y = R so that the radius of curvature is r = R/tan~ or 
and the normal unit load acting on the helical wound tow is 





Here another scaling factor was identified, size and helical angle which may relate to the lamina 
moduiii. 


E-t, G 


12 


or/'S 


for 



(36) 
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And here again the physical evidence would be an increase in laminate thickness ever that predicted 
but the weight would not increase, voids would. 

There exists some evidence that might be attributed to the above phenomena. Seventeen quarter 
scale specimen (Tables 7 and 8) thicknesses were measured to be 4 percent less than calculated while ne 
full scale (Table 10) shell measured 10 percent thicker than predicted. This is not a sound statistical base 
for modeling but neither is there a unique method at this time for relating thickness differences to resin 
elastic degradation nor how to distribute it to gel* and normal load events and to helical and hoop laminae. 

One clue for reducing the problem was detected in a fragment from a full scale burst test. The hoop 
lamina was noted to have a uniformly varying content of resin, and while occasional tows could be identified 
in most of the acreage, no voids between them were seen. On the other hand, helical tows were densely 
coated with resin and the gaps between them were noted with few exceptions. Furthermore, helical tows 
were not completely bonded to the adjacent layers which made it exceedingly easy to peel off each tow at 
will. These conditions were reported to have degenerated with article size as suggested above. 


Another clue was the strain analysis of hoop and helical lamina for resin dependency. The hoop 
lamina principal strains are coincident with the laminate strains and therefore defined by the g| character- 
istic. For a hoop fiber ultimate strain [equation (7)1 of 1.4 percent, the worst case (gj = 4.77) transverse 


strain was about 0.3 percent and marginal (Table 3) 
second order axis transformation was used 


1 e ‘ 1 

e x | 

| e 2 f = [Tl 

e 0 1 

I 7,2 1 

7x0 | 

~R~ 



to yield the transverse strain 

ei = eg (m~ t-n^/gj) (37) 

and shear strain 

7l2 = 2e 0 mn(l - 1/gj) (38) 


To determine the helical lamina principal strains, a 


HELICAL 



What this analysis implied was that even for large helical angles approaching hoop, the lamina 
transverse ultimate strain is reached long before the hoop fiber strain capability is achieved (Fig. 7). It does 
not, however, mean that the transverse modulus is zero. The 4-in. lamina tube only demonstrated that the 
3-in. length specimen separated at the weakest link when a strain of 0.293 percent was reached but the 
remaining length had additional straining capabilities. Recalling how gaps and tow densities (Fig. 6) may 
vary in a lamin a, and particularly the helical lamina, the transverse strain and associated modulus may be 

*gel refers to consistency or vitrification. 
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assumed to act in a loose and variegated pattern across the layer. This is not crazing as many suspected 
because it is not a uniform separation of fiber from resin within and out the tow. 

The full significance of this discussion begins to penetrate when one realizes that here lies the dis- 
tinction between wound and cloth graphite-epoxy lamina; the wound transverse modulus degenerates to a 
fraction of the classical, primarily through manufacturing and then through operational strains. 

And what about the inplane shear behavior of equation (38). That strain, on the average, is real 
too but the modulus is variegated like the transverse one. But unlike the transverse modulus, the helical 
tows may be visualized as a lattice network reaping substantial rigidity from inter-lamina shear between 
helical and hoop layers. Accordingly, helical lamina shear model should be expected to relate to both 
g-charactcristics. 

The final clue was deduced from the sensitivity analysis which clearly indicated that resin connected 
modulii cf the hoop lamina had only third order effects on gj and g-» constraints as determined from the 

preliminary model. It has even less effects because of the disputed model so why bother revising E 2 , Gj 2 , 
and vj 2 models as presented in equations ( 1 7)? 

Having disposed of insignificant hoop lamina properties, all geometric properties of the six sets of 
specimen were listed. Table 16, that may be useful to distinguish them in the selection of scaling factors. 

TABLE 16. SPECIMEN GEOMETRIC PROPERTIES 


Parameter 

warn 

DV-36 

‘ " ’ - — 

ABA2 

DA01.2 

STA. Jr. 

DA(o 

Scale, S 

0.25 

0.25 

0.138 

1.0 

0.25 

1.0 

Helical, a 

0.506 

0.506 

0.506 

0.556 

0.585 

0.585 

Thickness, t 

0.301 

0.300 

0.181 

1.168 

0.314 

1.308 

Layers, N 

13 

13 

8 

43 

19 

49 

Hoop. H 

0.302 

0.303 

0.277 

0.229 



0.245 

0.221 


Once the environments affecting the transverse modulus of the helical lamina were explained, it was 
relatively easy to arrive at a reasonable reduction. If the hoop strain is the principal mechanism exercising 
the helical transverse modulus, then the spacing of helical tow remnants attached along the hoop tows should 
relate to the proportion of helical lamina engaged. Then the gaps and lack of resin coalescences between 
helical tows should account for the balance of reduction. Upon review of the fragment from full-scale burst 
test, tattle-tail epoxy spotting and helical tow scraps on both sides of the hoop lamina were visually 
estimated at 25 percent. After some manipulations and iterations of parameters, a reasonable prediction of 
transverse modulus of the helical lamina was empirically obtained, 


E? = 0.27 E 2 V(« 2 /NS) 


( 39 ) 


where E 2 is the transverse modulus of the hoop lamina and is calculated from equation (17). 
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Conceding that the transverse modulus had diminished, then the helical lamina major Poisson’s ratio 
must likewise be expected to degrade. In the absence of explicit evidence, and because of its relative insig- 
nificance, the helical lamina Poisson’s ratio was arbitrarily modeled. 


a 

"12 


= 0.53 Fj2 [« 2 /NS] 


0.04 


(40) 


The helical lamina shear modulus was perhaps the most interesting to modify. It had the property 
to skew the g-characteristics (reduce one and increase the other simultaneously) and establish a precept to 
dispersions noted in Table 1 4. The inplane shear compliance to g-characteristics and parameters in Table 1 6 
was rationalized as follows. The modulus decreases as g-> decreases and. to a lesser extent (not a continuous 

lamina), as gj increases, was concluded from the sensitivity study. The modulus increases when the inter- 
lamina surfaces increase, as discussed above, which may be related to the number of layers (N). However, 
when the specified layer thickness is small (t/N) and uncontrolled by inadequate compaction, then inter- 
lamina contact decreases and shear rigidity decreases. Modulus particularly is reduced with increased scale 
represented by the thickness. These parameters were cast such that resulting values would allow for reason- 
able adjustments in longitudinal modulii of hoop and helical lamina when combined into equations (11) 
and (21) to agree with hydrotest results of Table 12. The final inplane shear modulus of the helical lamina 
was expressed by 


g“ 2 = 0.0375 G, 2 ^- y/ (N/Hgj) (41) 


Having resolved six of the eight lamina elastic constants, there was little latitude left to manipulate 
the two fiber oriented modulii, hoop and helical. The hoop lamina modulus turned out. 



0.089 


Ejg2 

tH 



(42) 


was dominated by the g-> characteristic as expected from the sensitivity study and was weakened by the 
total hoop layer thickness. There was a hint of reduction with increase of gj characteristic and scale but it 

was suspected that this radical would vanish when more full-scale hydrotest and dog-bone data becomes 
available. A remarkable outcome of this hoop lamina expression was that the longitudinal modulus of the 
ABA2 and STA Jr. specimen were predicted greater than limit equation (17) and represented by Ej in 

equation (42). Both specimens were subscale : ABA2 was smallest and STA Jr. had the most layers of all 
quarter-scale articles. One explanation may be that the smallest articles allowed for better compaction and 
the STA Jr. allowed for more compaction opportunities, both leading to reducing the resin content of equa- 
tion (13). Compaction is applied after each hoop wind with a heavy cloth belt to minimize void. 

Finally, the longitudinal modulus of the helical lamina was formulated 
a ^1 ■> 

Ej = 0.24 — [S 3 N (1.308/t)] , (43) 


ORIGINAL PACZ ® 
OF POOR QUALITY 
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H 

as a function of the hoop lamina longitudinal modulus E j . Scaling with respect to diameter (S) tends to 

drop out should the radical reconcile to a cube root through a later update. The modulus was inversely 
responsive to scaling of the thickness (1 ,308/t). 

Curve fit of lamina elastic models with specimen sets are presented in Figure 8. which clearly demon- 
strates the scaling phenomena (correlation is better than 0.9). Since the slope was sensitive to full-scale 
articles and since that data is the least available at this time, a later update is incumbent and should 
incorporate results from scheduled dog-bone specimen tests of full-scale laminate. Curve fits decisively 
demonstrate scaling dependency of lamina elastic properties. Scaling factors might have been more purely 
defined if informal test data were not so parsimoniously released by the Prime Contractor. But more to the 
point is that scaling-up article size decreases lamina stiffness potential. It would seem that a weight decrease 
of 3 percent or more might be realized through a modest development of tow and resin deployment and 
control at the winding feed point. The scope is clear. 




vaui 

ADV-3S 
OABA2 
e DM1.2 
O ST A JR 
□ DA03 



Figure 8. Scaling dependency of lamina elastic constants. 
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Equations (17) were adjusted by scaling equations (39) through (43) to calculate the lamina con- 
stants and laminate stiffness for all the specimen sets and are listed in Table 17. The resulting laminate 
stiffnesses represent an adjusted set of estimated parameters that were recycled into equations (25) and 
(27) to obtain a more fitting set of “observed” parameters (Table 18) and as manifested by negligible errors 
listed in Table 19. Compare these with Table 14. 

TABLE 17. ADJUSTED ESTIMATES OF LAMINA AND LAMINATE PROPERTIES 


Parameters 

QS1.2 

DV-36 

ABA2 

DA 1, 2 

STA. Jr. 

DA03 

Helical Lamina 







E“ Msi 

16.6 

17.3 

20.1 

17.2 

19.7 

16.9 

E“ Msi 


0.16 

0.20 

0.10 

0.15 

0.10 

Msi 

0.51 

0.52 

0.41 

0.44 

0.30 

0.42 

„« 

12 

0.14 

0.14 

0.15 

0.13 

0.14 

0.13 

Hoop Lamina 







E 1 * Msi 

18 

18.4 

22.4 

18.2 

20.1 

18 

E 2 Msi 

1 

1 

1 

1 

1 

1 

G j 2 Msi 

0.7 

0.7 

0.7 

0.7 

0.7 

0.7 

"12 

0.3 

0.3 

0.3 

0.3 

0.3 

0.3 

Ajj Matrix 
Aji x 10 6 

2.157 

2.226 

1.63 

8.709 

2.337 

8.994 

A 22 x 10 6 ' 

1.902 

1.937 

1.329 

6.540 

2.022 

7.413 

A 12 x 10 6 

0.560 

0.579 

0.456 

2.860 

0.913 

3.349 

A 66 x10 6 

0.693 

0.714 

0.526 

3.348 

1.007 

3.876 

Dy Matrix 

x 10 3 

x 10 3 

x 10 3 

10 6 

10 3 

10 6 

D 11 

17.66 

18.01 

5.28 

1.014 

19.60 

1.289 

L 22 

10.24 

10.31 

2.05 

0.733 

13.49 

1.096 

*>12 

4.42 

4.53 

1.48 

0.321 

7.33 

0.462 

D 66 

_ 

5.38 

5.49 

1.67 

0.376 

8.06 

. 0.536 
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TABLE 18. ADJUSTED OBSERVED LAMINA PROPERTIES 


Parameters 

Q1.2 

DV36 

ABA2 

DAI. 2 

STA Jr. 

DA03 

A 11 

2.153 

2.226 

1.630 

8.804 

2.336 

8.992 

A 22 

1.883 

1.943 

1.325 

6.565 

2.005 

7.362 

a 12 

0.548 

0.581 

0.455 

2.871 

0.905 

3.314 


TABLE 19. DEPARTURE OF SCALED “ESTIMATED” FROM “OBSERVED” PROPERTIES 

(Percent) 


Parameters 

Ql,2 

DV36 

ABA2 

DAI ,2 

STA Jr. 

DA03 

A 11 

0.09 

0 

0 

-1.08 

0.04 

0.02 

a 22 

0.90 

-0.31 

0.30 

-0.38 

0.85 

0.69 

A 12 

2.19 

-0.34 

0.22 

-0.38 

0.88 

1.05 


IX. DYNAMICS RESPONSE DILEMMA 


Two quarter scale specimen of the type investigated here were dynamically tested and compared 
with modal predictions using unsealed stiffness matrices. Both shells were vibrated to experimentally deter- 
mine axial, torsional and shell frequencies. Axial and torsional test results compared favorably with predic- 
tions. However, DV and STA Jr. shell frequencies were 1 5 and 8 percent lower, respectively, than predicted 
which implied a weakened Dj 2 stiffness component as well as a scaling consideration. Differences in 

laminate bending have been amply published [10,111 and will be investigated for the FWC. However, the 
lamina properties disparity cause and scaling for static conditions was paramount. 


X. CONCLUSIONS 


It was inevitable that high performance graphite-epoxy would be counted among space-age materials. 
Many unexpected singularities have been identified and mastered during the FWC development, such is the 
wonder of this material, such was the experience with laminate stiffness phenomena. 

Graphite-epoxy wound laminate stiffness is unlike its cloth laminate kin in that the epoxy con- 
nected lamina properties cannot be fully realized. The mandrel curvature and gel condition of the resin are 
primary causes affecting the coalescences of tow to adjacent tows and layers to adjacent layers. These condi- 
tions were defined as scaling factors related to distinct geometric parameters, though polymer chemists might 
select more appropriate time factors: they should. 
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Scaling was amply supported by test evidence, reported observations and curve fits. Having disclosed 
the phenomena and linking it to stiffness reduction should provide a basis and purpose for winding improve- 
ments to realize about three percent weight savings and more efficient structural performance. Ten percent 
reduction is possible. 

In the meantime, serious effort should be directed to translating the “A-basis” ultimate stress 
(strength) of laminates from quarter-scale specimen to full-scale flight articles with the new baselined system. 
Results from this investigation provide the insights and rational basis for it. 
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APPENDIX A 


SENSITIVITY ANALYSIS 


Sensitivity analysis is a formidable tool in aerospace design to identify and concentrate on the more 
significant performance variables. Because of the multitude of independent variables that must be satisfied 
in the design of a high performance filament wound pressure bottle, the author suggested the technique 
during the FWC feasibility study and imposed it on this contracted project to insure that: 

1) Significant independent variables affecting stiffness were identified 

2) Stiffness performance dispersions were predicted. 

3) A criteria basis was established to specify material acceptance and manufacturing tolerance 
controls. 

During this study, sensitivity analysis again proved its versatility as a diagnostic method for gaining 
insight into lamina in-situ properties. 

Lamin .;e membrane characteristics of interest were defined by 


_ 2A 11 ~ a 12 
§I a 22 ~ 2a 12 

_ A 11 a 22 ~ A 12 
2 a 11' A 12 


( 21 ) 


which provide the starting noint of this analysis. The variables of concern are the orthotropic elastic 
properties of the hoop lamina and the helical lamina which are assumed to be different. Using the chain rule 
of differentiation on the two characteristic equations, the sensitivity coefficient derivations follow. 

A) Sensitivity of g-characteristics with respect to laminate properties: 


dgj - b j 3A 1 1 + b2 9 -^22 * ^ Aj 2 


dg 2 = b 4 3 Aj j + b 5 d A 27 + bg 3 Aj 2 


(A-l) 


where 


b-7 = 1 /(A22 " 2A 12^ 


b 8 = l/(2Aj j - A 12 ) 


(A-2) 
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and 
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bj = 2b? , b^ - gjb7 , b3~(2gj-l)b7 

A-3) 

b 4 ~ * A 22 ~ 2g 2* b 8 * b 5 = A Il b 8 ’ b 6 = ^2 _2A 12^ b 8 
B) Laminate Stiffness with respect to lamina variables: 

Using expiession (9), the extensional stiffness properties of the laminate may be simplistically formu- 
lated from equations (8) and (10) to provide a convenient set of lamina related properties for this analysis. 

Assuming — I. and designating hoop properties with a bar, then 

Aj j = t(l~H) (Ejm 4 + E2(2vj2m^n^ + n 4 ) + 4Gj2m 2 n 2 l + (t H E-i) 

A 22 = t(l-H) (Fjn^ + E2(2i»j2m^n^ + m 4 ) + 4G j 2m 2 n 2 1 +(t H Ej) (A-4) 

A l 2 = t(l-H) (Ejm-n- + E-»(m 2 n 2 ♦ r'pOn 4 + n 4 ))-4Gj2m 2 n 2 ] + (t H Pp E2) 
and 

dAj j = cj 3Ej + c-» 3E-* + C3 3vj 2 + C4 3Gp + c$ 3E2 

dA">'> = dj 3Ej + d2 3E-> + C3 3vj2 + C4 3Gp + c 3 3Ej (A-5) 

dAp = f j 3E j + f2 3E2 + (3 dvp + (4 3G p + fj 3E2 + fg 3up 

where 

cj=t(l-H)m 4 , C2 = t(l-H)(2»'|2m^n“ + n 4 ) 

c 3 = t(l-H)2E2m 2 n 2 , c 4 = t(l-HM4m 2 n 2 ) , C5 = tH 

dj = t(l-H)n 4 , d 2 = t(l-HK2»' 1 2m 2 n 2 + m 4 ) 

(A-6) 

fj = t(l-H)m 2 n 2 , f2 = t(l-HXm 2 n 2 + t'i2(m 4 + n 4 )) 

f 3 = t(l-H)C 2 (m 4 + n 4 ) , f 4 = -c 4 , f 5 = tHv I2 , f 6 = tHE 2 
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OF POOR QUAL; f i 


C) The desired sensitivity coefficients are: 


dg|/dEj = bjCj + b 3 di + b_jf| 
dgj/dE-> = bjC-* + bjdj + b 3 f 3 
dgj/d»'l2 = bjc 3 + b 2 d 3 + b 3 f 3 
dgj/dGj2 = c 4 ^ b l + b 2 “ b 3^ 
dgj/dEj = b 3 c 5 
dgj/3F2 = bj c 3 + b 3 f 3 
dgi/dvj2 = b 3 f 6 


3g2/3Ej - b^j + bjdj + b 6 fj 
dg2/3E-» = b^C2 + bjdi + b^f-* 

Hji^M = b 4 c 3 + b 5 d 3 + b 6 f 3 
dg2/3G j2 = C4 (b4 + b 3 - bg) 
dg2/3E | = b 3 cj 
dgj/^2 ~ b 4 c 5 + b 6 *5 

dg 2 lto\2 = b 6 % 


Results for 10 percent change in lamina variables are given by 


characteristics. A% Change = 10% Ej(d^/3Ej)/gj 


(A-7) 


laminate 


, A 7c Change = 10% EjOAy/dEjVAy 


(A-8) 



rr *• • ... • 

jF POOR QUALITY. 

rfm OS 
Run 

INPUT H£L ANC . THICKNESS. H-RAT 10 


29. 3016 

. 3025 







XC FOR 10X 

VARIABLE 

VARIABLE 

NOn VALUE 

Cl COEFF 

C2 COEFF. Cl CHANCE 

C2 CHANCE 

MElICAL 

LAMINA 




Ei 

1 860E+00I 

8 33 IE-001 

9 766E-004 3 034E+00I 

1 588E-001 

E2 

1 000E+000 

2 268E+000 

7 800E-002 3 8G4E+000 

6 907E-001 

Vi 2 

3 000E-001 

2 S25E+000 

3366E-002 1 377E+000 

8 827E-002 

Cl 2 

7 000E-00I 

-5 429E-001 

9 278E-002 -G S42E-00I 

5 G78E-00! 

HOOP LAMINA 




El 

1 . 860E+0C > 

G . 883E-00I 

5 30 IE -002 2 238E+00I 

8 620E+000 

E2 

1 000E+000 

G 175E-001 

-4 24 IE -003 I .079E+000 

-3 707E-002 

vi 2 

3 000E-001 

1 . 256E+000 

-5 521E-003 6 587E-00! 

-1 448E-002 

LAMINATE 





A! 1 

2 52! E +000 

2 G37E+000 

-2 833E-002 1 I 62E+001 

-6 243E-00I 

A22 

2 1 65E+000 

7 544E+000 

5 8I0E-001 2 855E+001 

1 . 1 00E+001 

A12 

7 032E-00! 

1 377E+00! 

-G 05 IE-002 l G92E+00I 

-3 720E-001 

XCHANCE 

IN Aij per 10X 

change IN E 

• J 


VARIABLE 

AI • 

A22 

A 1 2 


HEL I CAl 

LAMINA 




El 

9 082E+000 

9.985E-001 

1 000E+00I 


€2 

5 36!E-0ei 

S 735E-001 

1 1 1 3E+000 


vi 2 

9.002E-002 

1 048E-001 

5.747E-001 


02 

4 201E-00! 

4 892E-001 

-1 586E+000 


HOOP LAMINA 




El 

0 000E+000 

7 839E+000 

0 000E+080 


E2 

3 618E-001 

0 000E+000 

3 892E-00I 


VI2 

0 000E+000 

0 000E+000 

1 . 1G8E-001 



RFM OV 
RUN 

INPUT HEL. ANC. .THICKNESS. H-RAT I C 
29. 3001 . 3029 

XC FOR 10X VARIABLE 

VARIABLE NOH. VALUE Cl COEFF . G2 COEFF. Cl CHANCE C2 CHANCE 
HELICAL LAOINA 

El I . 8S0E+001 9 267E-001 9.726E -004 3 025E+001 1588E-001 

E2 I 5500E+000 2 249E+000 7 859E-002 3 84 7E +000 6 898E-001 

VI2 3 000E-001 2 604E+000 3 352E-B02 I 37 IE +006 8.825E-002 

GI2 7.000E-001 -5 36 IE-001 8 218E-002 -6 586E-001 5 663E-001 

HOOP LAOINA 

El 1 8S0E+00! 6 840E-001 5 282E-002 2 233E+00I 8.622E+000 

E2 1 000E+000 6 145E-B01 -4 224E-003 1 078E+000 -3 708E-002 

V12 3 000E-001 l 248E+000 -5 464E-903 6 570E-001 -1 438E-002 

LAMINATE 

All 2 507E+000 2 641E+000 -2844E-002 1.162E+001 -6 2S8E-00! 

A 22 2 1 56E+008 7 525E+000 5 810E-80I 2 847E+001 I 088E+001 

A12 G 094E-001 1 373E+001 -G . 010E-002 I.68SE+801 -3.688E-00I 

XCHANCE IN Ai j per 10X change IN Ei j 

VARIABLE All A 22 AI 2 

HELICAL LAMINA 

El 8.081E+000 9 970E-00I 1.000E+00! 

E2 1 36 IE-001 6.725E-00I 1.I13E+000 

VI2 8 00 IE -002 1 . 047E-00I 5.747E-00I 

C12 4 200E-001 4 885E-001 -1.506E+000 

HOOP LAMINA 

El 0 000E+000 7 842E+000 0 000E+000 

E2 3 625E-00I 0.000E+000 3.898E-001 

V 1 2 0.000E+000 0 000E+000 1 . 1 70E-001 
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t'J 


OF 


A! JT/ 


REM ASA 
RUN 

INPUT HEL. ANG . . THICKNESS. H-RAT 10 
20. 181. 2773 

XG FOR 10X VARIABLE 


VARIABLE 

NOM VALUE 

Gl COEFF. 

G2 COEFF. Gl CHANGE 

G2 change 

HELICAL 

LAMINA 




£1 

1 860E+001 

t 53 IE +000 

5 547E-004 3 74OE+00I 

1 .GIBE-001 

E2 

i 00BE«eee 

4 01GE+000 

4 80GE -002 5. 28OE+000 

7 4Q8E-001 

v«2 

3 BflBE-eei 

4 63SE+000 

1.O16E-002 1 831 £+000 

8 . 970E-002 

GI2 

7 000E-00I 

-1 22 IE +000 

G 06OE-002 -I 126E+000 

6 G2QE-001 

HOOP LADINA 




El 

1 . 860E+00I 

1 077E+000 

2 91 6£ -002 2 638E+001 

8 . 463E+000 

E2 

1 000E+000 

8 873E-00I 

-2 . 358E-003 1 1 69E+000 

-3 G79E-002 

VI2 

3 000E-001 

2 012E+000 

-4 20 IE-003 7 O50E-001 

-2 00OE-B02 

L AH I NATE 





A 1 1 

1 561E+000 

5 652E+000 

-2 1 34E-002 1 162E+001 

-5 1O7E-001 

A22 

1 224E+000 

2 J 46E»00I 

5 810E-001 3 4G0E+001 

1 . 1 10E+001 

AI2 

4 353E-001 

4 009E+001 

-8 549E-002 2 298E+001 

-5 807E-001 

xchange 

IN A i j per 1 0X 

chcnge IN E 

' J 


variable 

A 1 1 

A22 

A 1 2 


HELICAL 

LAMINA 




El 

9 ! 20E *000 

1 098E+000 

1 005E+001 


E2 

1 367E-001 

7.404E-001 

III 8E+000 


V12 

9.039E-002 

1 . 1 52E-001 

5 773E-001 


C12 

4 218E-001 

5.378E-001 

-1 513E+000 


HOOP LAMINA 




El 

0 000E+000 

7 624E+000 

0.000E+000 


E2 

3 215E-001 

0 . 000E+000 

3 459E-001 


V12 

0 000E+000 

0 - 000E+000 

1 038E-001 



REM OA i 
RUN 

inpu t hel. ang . thickness. h- rat io 


3! 85.1 

168. 2288 


XC FOR 10X 

VARIABLE 

VARIABLE 

NOM VALUE 

G! COEFF. 

G2 COEFF. Cl CHANGE 

G2 CHANGE 

HELICAL 

LAMINA 




E! 

1 860E+00I 

1 209E+001 

5 734E-003 1 105E+002 

2 . 838E-001 

E2 

1 000E+000 

3 257E+001 

2 86 IE-001 1 600E+00) 

7 814E-001 

V ! 2 

3 000E-00I 

3 . 735E*001 

1 202E-001 5 . 505E+000 

9 . 596E-002 

Gl 2 

7.000E-00I 

-l 591E+001 

5 .551 E-001 -5 473E+000 

i 034E+000 

HOOP LAMINA 




El 

I 860E+001 

6 . 890E+000 

1 609E-00! G 296E+001 

7 965E+000 

E2 

1 000E+000 

4 709E+000 

-l G69E-002 2 314E+000 

-4 441E-002 

V12 

LAMINATE 

3 000E-001 

1 . 344E+001 

-4 669E-002 190 IE +000 

-3 727E-002 

A 1 1 

9 675E+000 

2 533E+000 

-1 004E-002 1 204E+001 

-2 584E-001 

A 22 

7 354E+000 

2 . 578E+00! 

G 021E-001 9 3I6E+00’ 

I 1 78E»001 

A12 

3 282E+000 

5 030E *00 1 

-! 747E-001 8.11 IE+001 

-1 526E+000 

XCHANGE 

IN Ai j per !0X 

CHonge IN E 

' j 


VARIABLE 

A 1 1 

A 22 

A12 


HELICAL 

LAMINA 




£1 

9 0'6E*000 

1 . 767E+000 

1 026E+001 


E2 

1 844E-00? 

7 853E-00! 

1 . 044E+000 


V I 2 

1 . 1 22E-001 

1 477E-001 

4 . 925E-001 


Gl 2 

5 . 238E-001 

S 891E-001 

-1 544E+000 


HOOP LAMINA 




El 

0 . 000E+000 

G 759E+000 

0 000E+000 


E2 

2 . 7S" > E-001 

0 . 000E+000 

2 443E-001 


VI 2 

0 000E+000 

0 . 000E+000 

7 328E-002 
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ORIGINAL FAQl. r.. 
OF POOR QUALITY 


REM STA 
RUN 

INPUT HEL ANC. .THICKNESS, H-RAT 10 
35 5. 314. .245 

XC FOR 10X 

VARIABLE NOM. VALUE Cl COEFF . C2 COEFF. Cl CHANCE 


HELICAL LAMINA 
El I 860E+001 


E 2 
V12 
Cl 2 

HOOP LAMINA 
El 
E2 
V 1 2 

LAMINATE 
A I I 
A22 
AI2 


I . 300E+000 
3 . B00E-00I 
7 800E-00! 

1 860E+001 

1 000E+000 
3 00 0E -001 

2 . 402E-000 

2 1 25E+008 
9 077E-001 


5 542E+000 
I 331E-00I 
’ 536E+0B1 
-6 200E+000 


2.200E-003 
7 I 50E-002 
3 492E-002 
I 57IE-00I 


8 20 IE *001 
I 059E+001 
3 . 665E+000 
-3 453E+000 

4.6IGE-OS1 
1 825 E +000 
1 430E+000 

1 .2 33E+001 
6 857E+00! 
5 . 624E+001 


311 9E+008 4 743E-002 

2 293E-000 -5 672E-003 

5 990E+000 -1 4I5E-002 

6 452E+000 -I 854E-002 

4.054E*001 6 . 165E-00I 

7 786E+001 -1 839E-00I 


XCHANCE IN A i j per 
VARIABLE All 
HELICAL LAMINA 


10X 


El 
E2 
V'l 2 
Cl 2 
HOOP 
E 1 
E2 
VI2 


LAMINA 


877E+000 
1 7 IE-001 
255E-001 
854E-00! 

000E+000 

203E-00I 

000E+000 


cHonge IN E i j 
A 22 

1 . 925E+000 
8 81 IE-001 
1 418E-001 
6 616E-001 

6 . 732E+000 

0 000E »000 
0 . 000E 1-000 


A 1 2 

1 . 029E+001 
1 005E1000 
4.51SE-001 
- 1 549Ei000 

0 000E+000 

2.542E-00! 
7 S27E-002 


REM OA 3 
RUN 

INPUT HEL ANC. .THICKNESS. H-RAT IQ 
33 . 55, I 308. 2214 

XC FOR 10X 


VARIABLE 

Non value 

Cl COEFF. 

C2 COEFF. Cl CHANCE 

HELICAL 

LAMINA 



El 

1 860E1001 

2 607E1001 

8 690E-003 1 735E+002 

E2 

1 000E1000 

6 G14E*001 

2 982E-00I 2 . 366E+001 

VI 2 

3 000E-00I 

7 59 1 £i001 

I.387E-001 8 . 1 47E+000 

Cl 2 

7 000E-00I - 

-3 627E+001 

7 093E-001 -9 082E»000 

HOOP LAMINA 



El 

1 860E+B01 

1 ,3G2Ei001 

I.787E-001 9 059£i00l 

E2 

1 300E+000 

8 997E1000 

-2 085E-002 3.2I8E+000 

v 1 2 

3 000E-001 

2 G74E+001 

-6 142E-002 2 870E>000 

LAMINATE 




A 1 1 

1 026E100! 

3 364 E +000 

-8 377E-003 1 . 234E+00I 

A22 

8 379Ei000 

4 702E+00I 

6.171E-001 I409E+002 

A12 

3 892Ei000 

9 235E+001 

-2 12IE-00I I 28GEi002 

XCHANCE 

IN Aij per 1 0X 

CMonge IN E 

• j 

VARIABLE 

All 

A 22 

A12 

HELICAL 

LAMINA 



El 

8 909E+000 

2 l 09E+000 

! 032E+001 

E2 

2 1 90E-001 

7 410E-001 

1 .007E+000 

vi 2 

1 264E-00I 

1 547E-001 

4 . 51 9E-001 

0 2 

5 898E-001 

7 . 2 1 9E-001 

-1 . 554E+000 

HOOP LAMINA 



El 

0.000E+000 

6.428E+000 

0 . 000E+000 

E2 

2 823E-001 

0.000E+000 

2 232E-00I 

V 1 2 

0 000E+000 

0. 000E+000 

6 . G96E-002 


VARIABLE 
C2 CHANCE 


3 725E-001 
S 507E-00I 
0 . 535E-002 

I . 001 E 1000 

8 . 028E+000 

5 I8IE-002 
3 863E-002 

4.052E-001 
1 . 1 92£i00t 
l 520Ei000 


VARIABLE 
C2 CHANCE 


3 . 795E-00I 
7.000E-00I 
9 767E-002 
1 . IS6E1000 

7 . 804 E 1000 

4 896E-002 
4 326E-002 

2 0I7E-001 
1 214E1001 
1 938E+000 
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APPENDIX B 


COMPUTER CODES AND RESULTS 


0 Laminate Stiffness Calculations 
• Newton Iterative Method 
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OKlvali'; AL ^ 

OF POOR QUALITY 


FIND 13 

OLD 

LIST 

iae REM LAMINAT STIFFNESS MATRICES, AiJ.Dij 
118 SET DEGREES 
120 PRINT 

130 PRINT "INPUT COMMAND < MATERIALS, LAYERS, OR OUTPUT): "J 
140 INPUT C* 

150 IF C**"NATERIALS“ THEN 180 

160 IF Ct* "LAYERS’ THEN 358 

170 IF Ct-’OUTPUT" THEN 490 

180 PRINT "INPUT NUMBER OF MATERIAL TYPES: "5 

190 INPUT M 

200 DELETE El . E2. E3. E4, E5. C 1 .C2. C3. C4 

210 DIM E1<M),E2<M*.E3*H - ,E4(M<.E5CM >,Cl>'M>.C2<H).C3<H),C4<rt> 

220 PRINT "INPUT HhTL PROP •. El • J • . E2> J) . G12f J > ,»U2< J> " 

230 INPUT J, XI « X2« X2. “4 

240 El < J>-X1 

250 E2CJ)*X2 

260 E3» J »=X3 

270 E4< J)=X4 

288 E 5 C J > = E 1 > J > •■ < E 1 •' J ' - E 2 < J ) * E 4 ‘ J)T2> 

296 CKJ>*Elf J)*E5' J> 

300 C2« J>*E2< J)*E5<J> 

310 C3< :J>*E3<J> 

320 C4( J>*E4( J>*E2< J>*E5< J) 

330 IF J<M THEN 238 
348 GO TO 120 

350 PRINT "INPUT « OF LAYERS: <0 FOR TAPE INPUT)"; 

360 INPUT L 

365 IF L-8 THEN 5000 

370 DELETE T1.B1.M1 

380 DIM T1<L),B1(L).M1<L> 

400 T0-0 

410 PRINT "INPUT LAYER «,MAT'L », THICKNESS, AHGLE" 

420 INPUT I,X1,X2.X3 
430 N1CD-X1 
440 T1<I)*X2 
450 Bid >*X3 
460 T8=T0-*Tl< I ) 

470 IF I<L THEN 420 

472 PRINT "STORE DATA ON TAPE 7 ‘YES OR NO)“ 

474 INPUT A# 

476 IF A*="Y£S” THEN 5110 
488 GO TO 120 

490 DELETE 62. B3,B4. 65.01.02,03,04. 66.05. 06 
492 DIM 82 •' L > , B3 < L > . 84 <: L > . B5 ' L > . B6 v L ) 

494 DIM 0 1 < L > . Q2 ‘ L ' « 03 < L > « 04 < L ' . 05 1 L ) * 06 > L > 

500 A 1=0 
510 A2=0 
520 A3=0 
530 A4=0 
540 A5=0 
556 A6=0 
560 Zl=-T0-2 
570 D1=0 
580 D2=0 
590 D3=0 
600 D4=0 
618 D5=0 
620 D6=0 

630 FOR 1-1 TO L 
640 N-HKI) 

650 B2<I>*C0S<eid))t2 

660 B3< I >*SIN<B1 <1 ) >12 

670 B4d )»82d )T2 

680 B5d)-B3<I)T2 

690 B6<I )*C0S(B1 <I)>tSIN<BKI>> 

700 Ql(I>»Cl<N)*B4d>+2*(C4<M>42*C3<N)>*B2<I)*B3<I>+C2(N)*B5d 

710 02< I >«C1 <N)»B5< I >*2*<C4<N)42*C3CN> )*e2< I >*B3< I >*C2<N>*B4< 1 




729 Q3< I>*(C1 (N>+C2(N)-4*C3(N) >tB2( I )*B3< I )+C4(N)*(84( I )+B3C I ) ) 

1 30 Q4(I)*(C1 <N)+C2(N>-2*(C3(N>+C4<N> >)tB2< I )XB3( I >+C3(N)X<B4( l)+B5< I 

740 05 < I >*(C1 <N)-(C4(N>+2*C3(N>>>*86( I )*B2( I ) 

750 05( I >*Q5< I > + <C4(N>-C2(N>+2*C3(N> >*B6<I >*B3(I ) 

760 06( l>*<Cl(NWC4(N>+2fC3(N>>>*B6(n*B3( I * 

770 B6< I >*06< I )+<C4(N>-C2(N>+2*C3(N> >*B6< n*B2( I > 

780 Al»Al+01(I>*Tl<n 

790 A2*A2+02(I>*T1( I> 

800 03*03+03(1 >*T 1<I> 

310 04*04+04 ( I >XT1 < I > 

820 85*85+05 < I >* T I ( I > 

838 A6*A6+66( I 'tTl ( I > 

948 22*21+T1(I) 

850 T3*22t3-Z1T3 

860 Dl=Dl+0i(I>*T3/3 

878 02*D2+02(n*T3 3 

880 D3=D3+Q3( I >*T3-'3 

890 D4=D4+04( I >*T3'3 

908 D5*D5+Q5( 1 >tT3-'3 

910 D6*D6+Q6< I >XT3'3 

920 21*22 

930 NEXT I 

940 PRINT ' 0 MATRIX p MATRIX” 

950 IMAGE 3E,1X,3E.1X,3E,2X.3E,1X«3E,1X,3E 
960 PRINT USING 950: Al, A3.A5.DI.D3.D5 
970 PRINT USING 550: A3,A2, A6.D3,D2,D6 
980 PRINT USING 950:A5.A6.A4,D5,D6,D4 
990 REM 

1000 PRINT "LAMINATE THICKNESS* "!T0 
1010 F1*A1XA2-A3t2 
1020 F2-Fl'<A2*Te> 

1030 F3*F1' , <A1»T0> 

1040 F4-A3'A2 
1050 FS-A4XT0 

1060 PRINT ” Ex Es Gxy Ux»* 

1070 IMAGE 3E.3X.3E,3X,3E,3X,3E 
1080 PRINT USING 1070: F2, F3, F5, F4 
2000 PRINT "COMPARE COMMAND (YES, NG>: "1 
2010 INPUT 0* 

2020 IF D*="VES“ THEN 2040 

2030 IF D4*"N0" THEN 120 

2040 PRINT "INPUT A1 1 , A22, A12. gl , g2" 

2050 INPUT Hb,H7,H8.P6,P7 
2090 H1*A1*A2-A3T2 
2100 H2=A2-2*A3 
2110 H3=2*A1-A3 
2120 H4=H3''H2 
2130 H5=H1/H3 
2200 Pl*(fil-H6mw0 H6 
2210 P2=(A2-H7>*100 -H7 
2220 P3=(A3-H8 '*100''H8 
2230 P4= ( H4-P6 > X 1 00 'P6 
2240 P5=>'H5~P7)«180-P7 

2250 PRIMT “PERCENT PREDICTED OUEP TEST" 

2260 IMAGE 10A.3E 
2270 PRINT USING 2260: "All PI 
2280 PRINT USING 2260: "A22" . P2 
2290 PRINT USING 2260: "A12" ,P3 
2300 PRINT USING 2260:"9l“,P4 
2310 PRINT USING 2260: "*2" , P5 
2320 GO TO 180 

4990 REM XX INPUT FROM TAPE XX 
5000 PRINT "ENTER FILE i* 

5010 INPUT F 
5020 FIND F 
5030 INPUT 033 :H» 

5040 PRINT H* 

5050 INPUT 933:L 
5060 DELETE T1.B1.M1 
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5170 PRINT 033:L,T1.B1,H1 
5180 GO TO 120 

8000 REN FILE 14 ; OSli-2 (LAYERS? 

8010 REN FILE 15 : DU36A (LAYERS > 

8020 REM FILE 16: A6A2 (LAYERS? 

8030 REH FILE 17: DA162 (LAYERS > 

8040 REH FILE 18: ST A JR (LAYERS' 

8050 REH FILE 19: DA3 (LAYERS' 


FIND 13 
OLD 

RUN 

INPUT COMHAND (MATERIALS, LAYERS, OR OUTPUT?: HATERIALS 

INPUT NUHBER OF HATERIAL TYPES: 3 

INPUT HAT'L PROP *.£1( J?,E2( J).G12> J),U12(J) 

1,18*6*1, *7, .3 

2. 18. 6.1. . 7. .3 

3. 6. 6. 2. . 76. .28 

INPUT COMMAND (HATERIALS. LAYERS, OP OUTPUT?: LAYERS 
INPUT « OF LAVERS: (0 FOP TAPE INPUT.?0 
ENTER FILE * 

14 

OSIi.2 

INPUT COMHAND (MATERIALS. LAYERS, OR OUTPUT?: OUTPUT 
A MATRIX D MATRIX 

2.453E+008 6.717E-001 0.000E+000 1.997E-002 5.290E-083 -2.690E-004 

6. 717E-001 2. 1 19E+000 0.000E+000 5.290E-003 1.175E-002 -9.253E-005 

0.000E4800 0.000E+800 7.876E-001 -2.698E-084 -9.253E-005 6.117E-803 

LAHINATE THICKNESS* 0.3016 

Ex Ey Gxy Uxy 

7.427E+000 6, 416E+080 2.612E+000 3.170E-801 

COMPARE COMMAND (YES, NO/.’ YES 
INPUT A11,A22,A12, 9 1, 9 2 

2.425. 1.8969.. 4916. 4. 77. 1.0 

PERCEHT PREDICTED OUER TEST 
All 1.146E+000 

A22 1 . 171E+081 

A12 3. 663E+801 

9 1 1.444E4001 

9 2 1.210E+001 
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ORIGINAL >9 

OF POOH Qw/.-.i • 

RUN 180 

INPUT NUMBER OF MATERIAL TYPES: 3 

INPUT MAT 'L PROP *,E1 ( J>,E2< J>,G12( J),U12< J) 

1,16«&,.16,.51,.14 
2. 18. 1 * ■ 7f 1 3 

3, 6.6. 2, . 76, #28 

A MATRIX 0 MATRIX 

2. 157E+008 5.599E-001 0.000E+000 1.766E-002 4.417E-003 -2.476E-004 

5.599E-001 1 . 902E+000 0.000E+000 4.417E-003 1.024E-002 -B.857E-005 

0#e80E+000 0. 080E+000 6. 933E-001 -2.476E-004 -8.837E-005 5.380E-083 

LAMINATE THICKNESS* 0.3016 

Ex Ey Gxy Uxy 

6.606E+000 5.823E*000 2.29SE+000 2.944E-00I 

PERCENT PREDICTED OUEP TEST 
All -1 . 104E+001 

A22 2.453E-001 

A12 1 . 389E+001 

9l 6. 816E-801 

*2 9.067E-001 

INPUT NUMBER OF MATERIAL TYPES: 


FIN013 

OLD 

RUN 

INPUT COMMAND (MATERIALS, LAYERS, OR OUTPUT): MATERIALS 

INPUT NUMBER OF MATERIAL TYPES: 3 

INPUT HAT 'L PROP *,E1(J),E2(J>,G12<J),U12< J> 

1,18.6,1, #7, #3 

2. 18.6.1. . 7. .3 

3.6. 6. 2. . 76.. 28 

INPUT COMMAND (MATERIALS. LAYERS. OR OUTPUT); LAYERS 
INPUT • OF LAYERS: «'0 FOR TAPE INPUT-0 
ENTER FILE • 

15 

DU36A 

INPUT COMMAND (MATERIALS, LAYERS, OR OUTPUT): OUTPUT 
A MATRIX D MATRIX 

2.439E+990 6.677E-001 0.000E+000 1.966E-002 5.20SE-003 -2.65BE-094 

6.677E-00 1 2.11 1E+000 0.006E+000 5.205E-803 1.159E-002 -9.141E-085 

0.000E+000 0.0O0E+000 7.831E-001 -2.658E-004 -9.141E-005 6.020E-003 

LAMINATE THICKNESS* 0.3001 

Ex Ey Gxy Uxy 

7.422E+000 6. 424E+000 2.609E+000 3. 164E-001 

COMPARE COMMAND (YES, NO): YES 
INPUT A11.A22.A12, 91*92 

2. 420. 1.930.. 5420. 4. 96. 1.03 

PERCENT PREDICTED OUER TEST 
All 7.672E-001 

A22 8.240E+000 

A12 2. 319E+001 

91 9.472E+000 

92 8.429E+090 

INPUT NUMBER* OF MATERIAL TYPES: 
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ORIGINAL P:--M 
OF POOR QUALITY 


RUN 180 

INPUT NUMBER OF MATERIAL TYPES: 2 

INPUT HAT'L PROP *,E1(J),E2< J>,G12<J>, U12CJ) 


IS. 


3.676E-004 

1.093E-902 

5.264E-001 


1.20.1.. 2.. 41. .15 

2.22.4.1. . 7. .3 
A MATRIX 

1 . 631E+000 4.567E-001 

4.567E-001 1 . 329E+000 

3.676E-004 1.893E-002 

LAMINATE THICKNESS* 0.131 
Ex Ey 

8. 142E+000 6.638E+000 

PERCENT PREDICTED OUER TEST 
Ail 3.927E+000 

A22 4.801E-001 

A12 -1.576E+000 

91 -4. 129E-001 

9 2 3.656E-001 

INPUT NUMBER OF MATERIAL TYPES: 


D MATRIX 
S. 285E-003 
1.478E-003 
4. 531E-007 


1.478E-003 

2.054E-003 

1.347E-005 


Gxy 

2.908E+000 


Uxy 

3.435E-0B1 


4. 531E-007 
1 • 347E-005 
1.667E-803 


F1ND13 

OLD 

PUN 

INPUT COMMAND < MATERIALS, LAYERS, OR OUTPUT): MATERIALS 

INPUT NUMBER OF MATERIAL TYPES: 2 

INPUT MAT'L PROP »,E1<J),E2< J>,G12< J),U12< J) 

1.18.6.1. . 7. .3 
2,18.6, 1 • . 7, . 3 

INPUT COMMAND (MATERIALS, LAYERS, OR OUTPUT): LAYERS 
INPUT * OF LAYERS: <0 FOR TAPE INPUT>0 
ENTER FILE • 

16 

ABA2 

INPUT COMMAND (MATERIALS. LAYERS. OR OUTPUT): OUTPUT 
A MATRIX D MATRIX 

1.568E+000 4.378E-801 3.671E-004 5.078E-003 I.4I5E-003 4.524E-007 

4.378E-001 1 .238E+909 8.929E-B83 1.413E-003 2.121E-B03 1.100E-085 

3.671E-004 8. 929E-003 3.099E-001 4.524E-007 1.100E-005 1.612E-003 

LAMINATE THICKNESS- 0.181 

Ex Ey Gxy Uxy 

7, 804E+980 6. 12 1E+099 2.817E+990 3.559E-081 

COMPARE COMMAND <YES,N0>: YES 
INPUT All, A22.A12, 91,92 

1 . 369. 1 . 323. . 464. 6. 77. . 696 

PERCENT PREDICTED OUER TEST 
All -4.497E-002 

A22 -7. 824E+000 

A12 -5. 645E+000 

9 1 1.247E+001 

m9 —7 

INPUT NUMBER'oF MATERIAL TYPES: 
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ORIG’:' ’ *"• ‘ 

OF F'OG.< Qc.V-.- - 

RUN I 86 

INPUT NUMBER OP MATERIAL TYPES: 3 

INPUT MAT'L PROP *,E1(J),E2<J),G12<J),U12( J> 

■Hill 

1.17.3. . 16. .52. .14 

2. 18.4.1. . 7. .3 

3.6. 6. 2. . 76. .28 

A MATRIX D MATRIX 

2.226E+000 5. 798E-881 0.000E+000 1.881E-002 4.526E-003 -2.552E-004 

5. 798E-001 1 . 937E+000 0.000E+080 4.326E-003 1.831E-082 -9. 103E-805 

0.000E+000 0. 000E+009 7.145E-881 -2.552E-004 -9.103E-803 5.489E-063 

LAMINATE THICKNESS- 0.3001 

Ex Ey Gxy Uxy 

6. 840E+-009 5. 953E+000 2.381E+000 2.993E-001 

PERCENT PREDICTED OUER TEST 
All -8.81 1E+000 

A22 -S.475E-001 

A12 6. 96< E+000 

91 3. 723E-001 

92 -2.996E-001 

INPUT NUMBER OF MATERIAL TYPES: 


FIND 13 

OLD 

RUN 

INPUT COMMAND (MATERIALS, LAYERS. OR OUTPUT): MATERIALS 

INPUT NUMBER OF MATERIAL TYPES: 3 

INPUT MAT'L PROP *,El< J),E2< J),G12< J>,U12< J> 

1.18.6.1. . 7. .3 
2, 18.6,1, *7, .3 

3. 17. 4. 1.6. . 71. .28 

INPUT COMMAND (MATERIALS. LAYERS, OR OUTPUT): LAYERS 
INPUT • OF LAVERS: <0 FOR TAPE INPUT)8 
ENTER FILE • 

19 

DA 003 JULY 19 1984 

INPUT COMMAND (MATERIALS, LAYERS, OR OUTPUT): OUTPUT 
A MATRIX D MATRIX 

1.030E+001 3. 848E+000 8.000E+000 1.468E+000 5.307E-001 -2.848E-002 

3. 848E+000 8.549E+000 0.800E+000 5.307E-001 1.253E+000 -1.345E-002 

0.000E+000 0.000E+000 4. 367E+000 -2.848E-002 -1.343E-002 6.039E-001 

LAMINATE THICKNESS- 1.30796 

Ex Ey Gxy Uxy 

6. 550E+000 3. 437E+000 3.338E+000 4.501E-001 

COMPARE COMMAND (YES, NO): YES 
INPUT All,A22,A12,9l'»2 
10.273,7.363,3.230,20,3. 764 
PERCENT PREDICTED OUER TEST 
All 2.316E-081 

A22 1.608E+001 

A12 1.840E»0bl 

91 -1 . 794E+000 

92 1.617E+881 

INPUT NUMBER OF MATERIAL TYPES: 
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original page SS 

OF POOR QUALITY 


RUN 188 

INPUT NUMBER OP MATERIAL TYPES: 3 

INPUT MAT'L PROP •, El ( J>, E2< J) , G12< J>, U12< J> 

liU.9i.lf .42, .13 

2. 18.1. . 7. .3 

3. 17.4. 1.6. . 71. .28 

A MATRIX 0 MATRIX 

8. 994E+800 3.349E+000 0.000E*e00 1.289E+808 4.623E-881 

3. 349E+008 7.413E+888 0.000E+000 4.623E-881 1.096E+888 

0.000E+000 0. 0B0E+900 3. 876E+000 -2.698E-802 -1.286E-882 

LAMINATE THICKNESS* 1.38796 

Ex Ey Gxy Uxy 

5. 719E+888 4. 714E+000 2.963E+00G 4.518E-681 

PERCENT PREDICTED OUER TEST 
All -1 . 245E+801 

A22 6.577E-801 

A12 3. 847E+000 

9l 2 . 313E+008 

l2 6.523E-881 

INPUT NUMBEr’oF MATERIAL TYPES: 


-2.698E-882 
-1.286E-882 
5. 366E-981 


PIND13 

OLD 

RUN 

INPUT COMMAND (MATERIALS, LAVERS, OR OUTPUT): MATERIALS 
INPUT NUMBER OF MATERIAL TYPES: 3 
INPUT MAT'L PROP *,E1<J>,E2(J>,G12(J>,U12<J> 
i , 18. 6, 1 < . 7, . 3 

2. 18. 6.1. . 7. .3 

3. 6. 6. 2. . 76. .28 

INPUT COMMAND (MATERIALS, LAVERS, OR OUTPUT): LAYERS 
INPUT » OF LAVERS: <0 FOR TAPE INPUT)8 
ENTER FILE t 
18 

STA JR 

INPUT COMMAND (MATERIALS, LAVERS, OR OUTPUT): OUTPUT 
A MATRIX D MATRIX 

2. 347E+000 8.663E-801 0.008E+800 1.969E-062 6.964E-003 -2.273E-004 

8.663E-001 2.087E+000 0.000E+000 6.964E-803 1.428E-002 -1.070E-004 

0.000E+000 0.000E+000 9.872E-001 -2.273E-004 -1.070E-004 7.980E-083 

LAMINATE THICKNESS- 0.314 

Ex Ey Gxy Uxy 

6. 331E+000 3. 628E+O00 3.144E+000 4.131E-081 

COMPARE COMMAND (YES, NO): YES 
IHPUT All, A22, A 12, g 1,92 

2. 349. 2. 005.. 9043. 19. 35. 1.026 

PERCENT PREDICTED OUER TEST 
All -6. 441E-002 

A22 4.078E+000 

A12 -4. 226E+888 

fl ~4. 414E+881 

m9 5 Mattasg 

INPUT NUMBER’ OF MATERIAL TYPES: 


48 




FIND 13 

OLD 

FUN 

INPUT COMMAND 'MATERIALS. LAYEPS. OR OUTPUT): MATERIALS 

INPUT NUMBER OF MATERIAL TYPES: 3 

INPUT MAT 'L PROP i, El < J) , E2< J> . G1 2' J > , U12< J> 

1 • 18. o, 1 . . 7. . 3 
2. 18. 6j 1 . . 7> .3 
3, 17. 4, 1.6, .71, .28 

INPUT COMMAND 'MATERIALS. LAYERS, OR OUTPUT): LAYERS 
INPUT i OF LAYERS: <0 FOR TAPE INPUT>0 
ENTER FILE # 

17 

DA 1*,2 

INPUT COMMAND (MATERIALS. LAYERS, OR OUTPUT): OUTPUT 
A MATRIX t> MATRIX 

9. 745E+000 3. 253E+000 0.000E+000 1.130E+000 3.650F-001 -4.667E-003 

3. 253F+000 7. 452E+000 0.000E+000 3.650E-001 8.339E-001 -1.961E-003 

0.000E+980 8 . 000E+000 3.716E+000 -4.667E-003 -1.961E-003 4. 170E-001 

LAMINATE THICKNESS* 1.16801 

Ex Ey Gxy Uxy 

7. 127E+800 5.450E+000 3.182E+000 4.366E-001 

COMPARE COMMAND (YES, NO): YES 
INPUT All,A22,A12,gl,92 
9.723,6.368,2.020, 17.91,3.363 
PERCENT PREDICTED OVER TEST 
All 2.224E-001 

A22 1 . 346E+001 

A12 1.537E+001 

9 1 -4. 1 15E+080 

92 1.361EV001 

INPUT NUMBER OF MATERIAL TYPES: 
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ORIGINAL PAGE f9 

OF POOR QUALITY 


RUN 180 

INPUT NUMBER OF MATERIAL TYPES: 3 

INPUT MAT L PROP i.El<J>,E2< J>.G12< J>,M12< J> 

1.17. 2.. 19. .44. .13 

2.19.2.1. . 7. .3 


3, 17.4, 1. b. 
A MATRIX 
8.709E+888 
2.860E+909 
8.000E+000 


‘1 • .28 


9. 008E+000 
0.000E*090 
3. 348E+0O0 


2. 860E+800 
6.S40E+000 
0.000E+000 

LAMINATE THICKNESS- 1.16891 
Ex Ey 

6. 385E+000 4. 795E*000 

PERCENT PREDICTED OUER TEST 
All -1 . 843E+061 

A22 -4.245E-001 

A12 1.411E+000 

91 -9.424E-981 

92 -3.670E-981 

INPUT NUMBER OF MATERIAL TYPES: 


D MATRIX 
1.814E+088 
3.210E-801 


3. 219E-081 
7. 328E-991 


-4.493E-803 -1.916E-883 


Uxy 

2.966E+088 


Uxy 

4.373E-881 


-4.493E-083 
-1 . 916E-883 
3. 7S8E-881 


FIND 6 

OLD 

LIST 

7 FIND 6 

198 REM LAMINATE HYDROTEST PROPERTIES, NEUTON ITERATIVE METHOD 
118 REN PREDICTED: Bl -A 1 1 , B2-A22,63*A12. TEST : B5-9 1 , 86*92 
128 REM DETERMINE: Al-Al It .A2-A22 * ,A3=A12* 

298 PRINT ’INPUT PRED.A11.A22.A12. TEST 9 I, 92 ’ 

218 INPUT B1,B2,B3,B4,BS 

230 PRINT ’INPUT TRIAL all,a22,al2’ 

240 INPUT B7.B8.B9 

258 PRINT ’INPUT NUMBER OF ITERATIONS ’ 

251 INPUT K 
268 1-0 
389 LI -2 
310 Ml— B4 
328 N1*2*B4-1 

338 Cl— 2*B7+B4*B8-<2*B4-1>*B9 
358 L2-2SB5-B8 
368 M2— B7 
370 N2-2*B9-B5 

380 C2— 2*B5*B7»B7*B8«85*B9-B9t2 

488 L3— 2»B7*<2*B4-mBl*<2*B4-l>-2*B3-B4t35+B9*<2*B4*2> 

418 H3-B9*<B4+4>-B4*B3-B2*<2*B4-l >-4*B4*B5+2*B8*<2*B4-l> 

420 N3— 2*B4*B 1 -2*B4*B5-4*B2+B8* < B4*4 > «07*2* < B4+ 1 > 

438 C3“<2*B4-1 >*B7t2-B7*<Bl*<2*B4-l >-2*83-B4*B5>-2*<B4*l >*B7*B9 

448 C3-C34B9* < 2*B4*B 1 «2*B4*B5«4*B2 > -B9*B8* <B4M ) ♦B8* < B4*B3*B2* < 2*84- 1 > > 

458 C3-C3*B8*4*B4*B5- < 2*B4- 1 > *38t2-B4*B5* < B 1 +2 *B3+4*B2 ) 

588 D1>L1*<H2*H3-N2*M3>-L2*<H1*N3-N1*N3>*L3*<H1*N2-N1*M2> 

518 D2«C1*<N2*N3-N2*H3>-C2*<N1*N3-N1*N3>+C3*<H1*N2-N1*M2> 

528 D3«L1*<C2*N3-N2*C3>-L2*<C1*N3-N1*C3>+L3*<C1*N2-N1*C2> 

538 D4-tl*<H2*C3-C2*H3>-L2*<Nl*C3-Cl*H3>+L3*<Hl*C2-Cl*M2> 

548 X-D2/D1 
558 Y-03/D1 
568 Z-04/D1 

578 PRINT ’X.Y,2. ’»X,Y,Z 
608 B7-B7+X 
618 B8-B8+Y 
'28 B9-B9+Z 

638 PRINT ’TEST All, A22. A12 "jB7.B8.B9 
648 I-I+l 

658 IF KK THEN 390 
668 END 
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